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In  the  first  half  of  the  dissertation,  an  algebraic  approach  is 
developed  based  on  "augmentation"  for  linear  discrete-time  time-varying 
systems.  Using  this  approach,  it  is  shown  how  to  design  feedback  control 
systems  for  a broad  class  of  time-varying  systems  which  result  in  a 
closed-loop  system  matrix  that  is  equivalent  to  a constant  matrix  with 
arbitrary  assignable  eigenvalues.  The  theory  yields  controllers  having  a 
much  lower  order  in  general  than  that  obtained  in  past  work.  In  the 
second  half  of  the  dissertation,  a discrete-time  time-varying  version  of 
Sylvester's  resultant  matrix  is  defined.  This  definition  is  then  used  to 
check  coprimeness  of  polynomials  and  polynomial  matrices  with  time- 
varying  coefficients,  to  examine  reachability  of  a given  system,  and  to 
design  compensators  which  result  in  a type  of  assignability. 


CHAPTER  ONE 
INTRODUCTION 


Time-varying  systems  can  be  found  in  many  practical  applications, 
such  as  linearized  models  of  nonlinear  systems  about  a time-varying 
nominal  trajectory,  adaptive  aircraft  systems  with  time-varying 
parameters,  integrated  circuits  with  time-varying  components,  and  so 
on.  Thus,  there  is  a good  deal  of  motivation  for  developing  the  theory 
of  linear  time- varying  systems. 

The  first  part  of  this  dissertation  is  concerned  with  developing 
the  algebraic  theory  of  linear  time-varying  systems.  Our  emphasis  is  on 
feedback  control  system  design  that  is  applicable  to  a broad  class  of 
time-varying  systems. 

A part  of  the  past  work  on  feedback  control  of  linear  time-varying 
systems  is  based  on  the  minimization  of  a quadratic  cost  and  is 
specified  in  terms  of  Riccati  difference  equations.  For  example,  CHENG 
[1979]  and  KWON  and  PEARSON  [1978]  have  obtained  stabilizing  memory less 
feedback  laws  for  uniformly  reachable  time-varying  systems.  HAGER  and 
HOROWITZ  [1976]  and  ANDERSON  and  MOORE  [1981]  have  defined  the  notion  of 
stabilizability  for  discrete-time,  time-varying  systems  and  have  derived 
the  optimal  control  in  memory less  feedback  form. 

Other  approaches  to  linear  time-varying  systems  are  based  on 
algebraic  formulations  which  yield  systematic  methods  for  computing 
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feedback  control  laws.  For  example,  WOLOVICH  [1968]  has  obtained  a 
stabilizing  feedback  law  for  continuous  time-varying  systems  possessing 
a fixed  lexicographic  basis.  MORSE  and  SILVERMAN  [1973]  have  derived 
coefficient-assignability  results  on  the  control  of  index-invariant 
continuous  time-varying  systems.  These  previous  results  can  be  easily 
extended  to  the  discrete-time  case.  Using  the  concept  of  n-cyclicity, 
KAMEN  and  HAFEZ  [1979]  have  developed  an  algebraic  theory  for  discrete- 
time time-varying  systems. 

In  general,  approaches  based  on  algebraic  formulations  provide  more 
flexibility  and  insight  than  the  one  based  on  the  minimization  of  a 
quadratic  cost.  But  all  the  past  work  based  on  the  algebraic  approaches 
mentioned  above  treats  very  restricted  classes  of  linear  time-varying 
systems,  such  as  index-invariant  or  n-cyclic  systems.  Therefore,  it  is 
desirable  to  develop  an  algebraic  theory  of  feedback  control  for  a much 
broader  class  of  linear  time-varying  systems. 

In  time-varying  systems,  it  is  possible  that  the  system  is 
reachable  at  all  times  in  N steps  where  N is  strictly  greater  than  the 
system  dimension  n.  This  is  due  to  the  lack  of  a "Cayley-Hamilton"  type 
theorem.  Therefore,  it  is  not  difficult  to  see  that  N-step  reachability 
is  a much  weaker  condition  than  n-step  reachability  which  in  turn  is 
weaker  than  index-invariance  and  the  existence  of  a fixed  lexicographic 
basis.  Moreover,  N-step  reachability  is  weaker  than  the  notion  of 
uniform  reachability  (see  KWAKERNAAK  and  SIVAN  [1972]). 

KAMEN  et  al . [1985]  and  POOLLA  [1984]  have  developed  an  algebraic 
theory  for  N-step  reachable  time-varying  systems  using  a transfer- 
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function  approach.  In  this  previous  work,  the  authors  have  successfully 
generalized  polynomial  matrix-f raction  methods  of  the  linear  time- 
invariant  system  theory  to  the  linear  time-varying  case,  and  have 
obtained  coefficient-assignability  results  for  N-step  reachable 
discrete-time  time-varying  systems  by  dynamic  state  feedback.  Time- 
varying  polynomial  matrix-fraction  methods  used  in  this  work  are  based 
on  a framework  specified  in  terms  of  noncommutati ve  (skew)  rings  of 
polynomials,  formal  series,  and  formal  Laurent  series,  all  with 
coefficients  in  a ring  of  time  functions. 

In  this  dissertation,  we  present  a constructive  method  for  finding 
a dynamic  state  feedback  controller  when  the  linear  time-varying  system 
IS  N-step  reachable  at  all  times.  Although  KAMEN  et  al . [1985]  have 
obtained  a dynamic  state-feedback  controller,  we  believe  that  our  scheme 
for  finding  a controller  is  much  more  flexible.  Moreover,  our  procedure 
yields  a controller  having  a much  lower  order  in  general  than  the 
controller  derived  by  KAMEN  et  al.  [1985]. 

The  second  part  of  this  dissertation  is  concerned  with  developing 
the  theory  of  the  resultant  of  linear  time-varying  systems. 

The  resultant  is  a classical  method  for  determining  the  existence 
of  a common  factor  of  two  time-invariant  (i.e.  constant-coefficient) 
polynomials.  From  the  resultant,  an  algorithm  for  determining  a 
greatest  common  divisor  (gcd)  of  two  polynomials  has  been  derived  (see 
BARNETT  [1983]).  There  are  several  definitions  for  time-invariant 
polynomials;  for  example,  Sylvester's  resultant,  Macduffee's  resultant, 
and  Bezout's  resultant.  For  the  definitions  of  these  resultants,  the 
reader  may  refer  to  KAILATH  [1980]  and  BARNETT  [1983]. 
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The  resultant  was  first  generalized  to  the  case  of  time-invariant 
polynomial  matrices  by  ROWE  [1972].  Later,  ANDERSON  and  JURY  [1976]  and 
BITMEAD  et  al . [1978]  have  modified  ROWE's  resultant  and  obtained  an 
algorithm  to  find  a gcd  of  polynomial  matrices. 

The  resultant  has  been  used  in  system  theory  to  obtain  minimal 
state-space  realization  of  transfer-function  matrices,  to  study 
decoupling  zeros  and  uncontrollable  and  unobservable  modes  of  given 
systems,  to  obtain  the  pole-zero  structure  of  given  multivariable 
systems,  and  to  design  dynamic  controllers  (see  CHEN  [1984],  BARNETT 
[1983],  and  KAILATH  [1980]). 

In  the  time-varying  case,  there  is  a notion  of  coprimeness  of  two 
time-varying  polynomials  and  polynomial  matrices.  This  notion  is 
defined  in  terms  of  a time-varying  version  of  the  Bezout  identity  (see 
Ilchmann  et  al . [1984]).  But,  this  definition  of  coprimeness  is  not 
related  to  a notion  of  the  resultant.  In  this  dissertation,  we  will 
define  the  discrete-time  time-varying  version  of  the  resultant  and  we 
will  find  the  connection  between  coprimeness  and  the  resultant.  This 
resultant  can  also  be  used  to  check  reachability  of  a given  time-varying 
system  and  to  design  feedback  control  systems. 

In  Chapter  Two,  we  summarize  basic  definitions  and  results  from 
KAMEN  and  HAFEZ  [1979],  KAMEN  et  al . [1985],  POOLLA  [1984],  and  so  on. 
Those  definitions  and  results  include  the  skew  polynomial  ring 
structure,  n-cyclic  systems,  N-step  reachability,  control -canoni cal 
form,  and  so  on. 
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In  Chapter  Three,  we  present  an  algebraic  approach  based  on 
"augmentation"  for  constructing  a generalized  control-canonical  form  for 
N-step  reachable  time-varying  systems.  If  a given  discrete-time  time- 
varying  system  is  N-step  reachable  at  all  times,  then  it  is  possible  to 
turn  the  given  system  into  a system  having  a fixed  basis  by  simply 
augmenting  the  system  state.  This  augmentation  is  always  possible  as  a 
result  of  the  special  properties  of  the  ring  of  time  functions.  From 
the  fixed  basis,  a control-canonical  form  for  the  augmented  system  is 
constructed  in  this  chapter,  which  we  call  the  generalized  control 
canonical  form.  Since  the  augmentation  forces  the  system  to  have  a 
fixed  basis,  the  augmented  system  is  index  invariant,  and  the  discrete- 
time  versions  of  the  results  of  WOLOVICH  [1968]  and  MORSE  and  SILVERMAN 
[1973]  can  then  be  applied  to  N-step  reachable  systems.  For  example, 
the  coefficients  of  the  closed-loop  characteristic  polynomials  of  N-step 
reachable  systems  can  be  assigned  arbitrarily. 

The  generalized  control -canonical  form  leads  to  the  design  of  a 
state  feedback  gain  matrix  for  the  augmented  system  which  assigns  the 
coefficients  of  the  closed-loop  system  matrix  arbitrarily.  Therefore, 
the  feedback  gain  matrix  plus  the  augmentation  comprises  the  controller 
for  the  N-step  reachable  system,  and  we  have  the  result  that  we  can 
construct  a closed-loop  system  whose  system  matrix  is  algebraically 
equivalent  to  a constant  matrix  with  arbitrary  assignable  eigenvalues. 

This  assignability  result,  which  was  first  obtained  by  KAMEN  et  al . 

[1985]  and  POOLLA  [1984],  is  discussed  in  Chapter  Four.  Also,  in  this 
chapter,  the  augmentation  approach  is  applied  to  design  of  output 
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feedback  controllers  when  the  system  is  given  by  an  input-output 
difference  equation.  We  note  that  a different  augmentation  approach  in 
the  time-invariant  case  was  considered  by  BRASCH  and  PEARSON  [1970]. 

The  augmentation  presented  in  this  dissertation  is  an  augmentation  of 
system  state,  whereas  the  one  considered  by  BRASCH  and  PEARSON  is  an 
augmentation  of  the  system  input  and  output. 

In  Chapter  Five,  cyclizability  of  N-step  reachable  time-varying 
systems  is  discussed.  HEYMANN  [1968]  has  developed  a direct  method  to 
find  a state-feedback  matrix  without  using  a control  canonical  form  for 
multivariable  time-invariant  systems.  Roughly  speaking,  the  system  is 
cyclizable  if  there  is  a feedback  gain  matrix  so  that  the  closed-loop 
system  is  reachable  through  one  of  the  inputs.  It  is  well  known  that  a 
time-invariant  system  is  cyclizable  if  and  only  if  it  is  reachable.  In 
the  time-varying  case,  however,  the  notion  of  cyclizability  has  not  been 
completely  investigated.  In  this  chapter,  we  shall  discuss 
cyclizability  of  the  augmented  system  for  N-step  reachable  time-varying 
systems.  This  result  is  then  applied  to  the  construction  of  feedback 
control  systems  for  N-step  reachable  time-varying  systems. 

Chapter  Six  treats  the  first  half  of  the  resultant  theory  for  time- 
varying  systems.  In  this  chapter,  we  define  a time-varying  resultant 
matrix  and  this  definition  is  used  to  check  reachability  of  a given 
time-varying  system.  This  chapter  is  divided  into  two  parts.  In  the 
first  part,  we  discuss  the  resultant  matrix  for  the  SISO  case,  and  in 
this  case,  the  resultant  matrix  resembles  Sylvester's  resultant.  For 
the  time-invariant  case,  the  reader  may  refer  to  BARNETT  [1983].  The 
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second  part  of  this  chapter  generalizes  the  definition  of  the  SISO 
resultant  matrix  to  the  MIMO  case.  The  reader  may  refer  to  ROWE  [1972], 
ANDERSON  and  JURY  [1976],  and  BITMEAD  et  al . [1978]  for  the  time- 
invariant  MIMO  case. 

The  second  half  of  the  resultant  theory  for  time-varying  systems, 
which  includes  the  relationship  between  coprimeness  and  the  resultant, 
design  of  compensators,  and  so  on,  is  covered  in  Chapter  Seven. 

Finally,  in  Chapter  Eight,  we  make  some  concluding  remarks  and 
discuss  some  open  problems  related  to  the  topics  covered  in  this 
dissertation. 


CHAPTER  TWO 
BACKGROUND 

The  purpose  of  this  chapter  is  to  give  various  preliminary 
definitions  and  results.  It  is  intended  to  help  the  reader  in 
understanding  the  material  covered  in  this  dissertation. 

With  Z - set  of  integers  and  R = field  of  real  numbers,  let  A 
denote  the  set  of  all  functions  defined  on  Z with  values  in  R.  With 
pointwise  addition  and  multiplication,  A is  a commutative  ring  with 
identity  1,  where  l(k)  = 1 for  all  k in  Z. 

Let  a be  the  right-shift  operator  defined  as 

a : ^ —►A 

a(k)  ' (aa) ( k )=a (k-1 ) , for  all  integers  k. 

Here,  k is  the  discrete-time  index  and  the  operator  a is  a ring 
automorphism  on  A.  The  inverse  of  a is  the  left-shift  operator  on 
A.  The  ring  A with  the  ring  automorphism  a is  called  a difference  ring 

For  any  positive  integer  n,  let  R^  denote  the  vector  space 
consisting  of  all  n-element  column  vectors  over  the  field  R.  Then  in 

terms  of  the  above  constructions,  we  have  the  following  notion  of  a 
system. 
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(2.1)  DE_FINITION.  An  m-input  p-output  n-dimensional  linear  time-varying 
discrete-time  system  over  A is  a triple  (F,G,H)  of  matrices  over  A where 

F is  n X n , G is  n X m , and  H is  p X n , together  with  the  dynamical 
equations 

x(k+l)  = F(k)x(k)  + G(k)u(k) 
y(k)  = H(k)x(k) 

Here.  x(k)  e R"  is  the  state  at  time  k,  u(k)  e is  the  input  at  time 
k,  and  y(k)  e RP  is  the  output  at  time  k. 

For  any  positive  integer  n,  let  A^  denote  the  free  A-module 
consisting  of  all  n-element  column  vectors  over  A.  Given  an  n x n 
matrix  F over  A,  let  Sp  denote  the  operator  on  A*^  defined  by 

Sp  : Afi  — ► 

(2.2) 

v(k)  I ► F(av)  = F(k)v(k-1),  for  all  integers  k. 

Note  that  Sp(av)  = F0(av)  = F(aa)(av)  = (aa)Sp(v).  The  operator  Sp  is 
said  to  be  a semi  linear  transformation  (KAMEN  and  HAFEZ  [1979]). 

The  system  (F.G)  over  A is  said  to  be  reachable  in  N steps  at  all 
limes  (or  N-step  reachable  at  all  times)  if  for  any  state  x belonging  to 
n-dimensional  Euclidean  space  R^  and  for  any  time  k , there  is  a control 
sequence  (u(k-N) ,u(k-N+l ) , . . . , u(k-l)}  that  drives  the  system  from 
x(k-N)  - 0 to  x(k)  = X.  A system  (F.G.H)  over  A is  said  to  be 
observable  in  N steps  at  all  times  if  and  only  if  the  system  (F',H‘)  is 
reachable  in  N steps  at  all  times,  where  the  prime  denotes 
transposition. 
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Let  C[^(k)  denote  the  n x mN  reachability  matrix  defined  by 

Cn(I<)  = [G(k)  (SpG)(k)  . . . (Sp^“^G)(k)] 

It  is  well  known  (WEISS  [1972])  that  the  system  (F,G)  is  N-step 
reachable  at  all  times  if  and  only  if 


rank  C[^(k)  = n , for  all  integers  k. 


(2.3)  example . Consider  the  system  (F.G)  over  A given  by  F(k)  0 1 

and  G(k)  = 


Lod(k)  1 


where  od(k)  = 0 if  k is  even  and  od(k)  = 1 othe 


rwi se. 


Then  C3(k) 


1 

od(k) 


1 

1 


which  has  rank  2 for  all 


integers  k.  Thus  the  system  is  reachable  in  3 steps  but  not  in  2 steps. 
Therefore,  in  discrete-time  time-varying  systems,  the  system  (F,G)  can 
be  reachable  in  N steps  at  all  times,  where  N > n.  This  situation  may 

occur  since  we  do  not  have  a "Cayley-Hamilton"  theorem  in  the  time- 
varying  case. 

Let  the  reachability  gramian  W(kQ,ki)  Integers  kj  > kg  be  the 
symmetric  nonnegative-definite  matrix  defined  by 


W(ko.ki)  = I ^(ki,k+l)G(k)G'(k)$’(ki,k+l) 

k=kg 
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where  the  prime  denotes  transposition  and  where  $(ki,k+l)  is  the 

transition  matrix  of  the  system  given  by  (see  KWAKERMAAK  and  SIVAN 
[1972]) 

*(k,ko)  = I F(k-l)F(k-2)  . . . F(ko)  , k > kg+l 

^ ^ * k=ko 

Then  it  can  be  shown  that  W(ko,ki)  = Cfj(ki-l)C^(ki-l) , where  N = kj  - 
kg.  Therefore,  the  system  (F,G)  is  N-step  reachable  at  all  times  if  and 
only  if  there  exists  a fixed  integer  N such  that  W(ko,ki)  is  nonsingular 
at  all  times  kg  e Z with  kj  = kg  + N. 

Let  C(k)  be  any  p X q matrix  over  A such  that 

rank  C(k)  < n for  all  k e Z . 

Then  there  exist  a q x q matrix  M(k)  over  A and  an  p x n matrix  W(k) 

over  A with  rank  M(k)  = q and  rank  W(k)  < n for  all  integers  k,  and  such 

that 


C(k)M(K)  = [W(k)  : 0]  . 


This  is  the  discrete-time  version  of  Dolezal's  theorem  (see  FERRER 
[1984]). 

Let  us  consider  column  vectors  vi,V2.  . . . , v^  belonging  to  A^ 
where  q is  a positive  integer.  We  say  that  Vq  is  an  A-1 inear 
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combination  of  vj,  V2,  • . '^q-i  if  there  exist  aj , a2 

aq_i  e A such  that  Vq  = a.  v- . 

i=l 

Again  let  (F,G)  be  an  m-input  N-step  reachable  system  with 
reachability  matrix  Cf^(k).  Assume  that  rank  G(K)  = m for  at  least  one 
value  of  k.  This  does  not  result  in  any  loss  of  generality.  In 
particular,  if  rank  G(k)  < m for  all  integers  k,  by  the  discrete-time 
version  of  Dolezal's  theorem,  there  is  an  m x m invertible  matrix  M(k) 
and  an  n X m matrix  G(k)  with  m < m such  that 

G(k)M(k)  = [G(k)  I 0] 

with  rank  G(k)  = m for  at  least  one  value  of  k.  Then,  we  can  work  with 
the  pair  (F,G),  since  if  we  let 

M(k)  = [Mi(k)  I M2(k)] 

where  Mi(k)  is  an  m x m matrix  and  M2(k)  is  an  m x (m-m),  then  the 
system  (F,G)  is  equivalent  to  (F,G)  via  the  input  transformation  M^fk). 

With  the  columns  of  G denoted  by  gi,  g2,  . . . , g^,  consider  the 
array  of  vectors  shown  below: 

9l  92 

Sf9i  Spg2 

• • 

• • 

• • 

S C N-1 

^F  9l  Sp"  "g2 


9ni 

^F9m 
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By  searching  the  rows  of  this  array  from  left  to  right,  we  select  those 
vectors  that  are  not  equal  to  an  A-1 inear  combination  of  previously 
selected  vectors.  Vectors  that  are  A-1 inear  combinations  of  previously 
selected  vectors  are  discarded.  In  this  selection  process,  it  is  easy 
to  see  that  if  Sp^g^-  is  discarded  for  some  integers  i,j  (where  1 < i < m 
and  1 < j < N-1),  then  Sp^g^-  will  also  be  discarded  for  j < A < N-1. 
Hence,  the  selection  process  will  result  in  an  array  of  the  form 


9l  92 

Spgi  Spg2 


9m 

^F9m 


There  are  no  gaps  in  this  array;  that  is,  if  SpJg^  is  in  the  array,  then 
Sp^'gi  is  in  the  array  for  0 < r < j.  Since  the  given  system  was 

assumed  to  be  N-step  reachable  at  all  times,  the  vectors  in  the  above 
array  must  generate  A*^. 

In  this  work,  we  shall  refer  to  the  integers  n^ , n2,  . . . , n„,  in 
the  above  array  as  the  generalized  controllability  indices  of  the 
system.  If  the  system  (F,G)  is  time  invariant,  the  n^-  are  the  usual 
controllability  indices  of  the  system.  Letting  q = n^  + n2  + . . . + 
we  have  that  N < q < mN.  If  N = n (i.e.,  the  system  is  n-step 
reachable),  and  q = n , we  shall  say  that  the  system  (F,G)  has  a fixed 
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_ba£i_s_.  Note  that  any  reachable  time-invariant  system  has  a fixed 
basis.  If  the  system  (F,G)  has  a fixed  basis,  a reodering  of  the 
elements  of  the  above  array  corresponds  to  a fixed  lexicographic  basis 
as  derived  in  the  work  of  WOLOVICH  [1968]. 

A system  (F,G)  is  said  to  be  index  invariant  (MORSE  and  SILVERMAN 
[1973])  if,  for  each  i = 1,  2,  . . . , n,  rank  C^-(k)  = q^-  = constant, 
for  all  integers  k,  and  = n where  C^(k)  = [G(k)  (SpG)(k)  . . . 
(Sf-^G)(k)]. 

Any  index-invariant  system  and  any  system  having  a fixed  basis  are 
n-step  reachable  systems.  But  the  converse  is  not  true.  Any  system 
with  a fixed  basis  is  index  invariant.  But  an  index-invariant  system 
does  not  always  have  a fixed  basis.  Therefore,  the  existence  of  a fixed 
basis  is  a very  strong  condition  in  the  time-varying  case. 


(2.4)  EXAMPLE.  Let  F(k)  = 


"1 

1 

O 

and  G(k)  = 

"k 

O' 

.0 

i_ 

1 

O 

1 

’k  0 k-1  O' 

.0  k 0 k-1 


Then  C2(k)  = 

which  has  rank  2 for  all  integers  k.  Thus  the  system  is 


reachable  in  2 steps,  but  (F,G)  is  not  index  invariant  and  does  not  have 
a fixed  basis.  j--j 


Then  C2(k)  = 

which  has  rank  2 for  all  integers  k.  Thus  the  system 


EXAMPLE.  Let  F(k)  = 

"0  O' 

and  G(k)  = 

1— H 
+ 

1 

.1  1_ 

-0  0 . 

k k+1  0 0 

0 0 k-1  k 
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is  reachable  1e  2 steps.  Since  rank  G(k)  = 1 for  all  k.  the  system  Is 
also  index  invariant.  But  the  system  does  not  have  a fixed  basis.  [] 

From  the  above  discussion,  we  see  that  N-step  reachable  systems  are 
the  largest  class  among  the  other  classes  of  systems  mentioned  above, 

including  n-step  reachable  systems,  index-invariant  systems  and  systems 
having  a fixed  basis. 

Let  A[z]  denote  the  set  of  all  finite  sums  of  the  form  I z'a.  in 
the  indeterminate  z with  coefficients  a^-e  A written  on  the  right.  With 
the  usual  addition  and  multiplication  defined  in  terms  of 

z^zJ  = zT+j  , i.j  e z.  i.j  >0 

(2.6) 

az  = z(aa)  , a e A , 

ACz]  is  a noncommutative  ring  called  a (right)  skew  polynomial  ring  or  a 
polynomial  difference  ring. 

A semilinear  transformation  Sf  defined  by  (2.2)  induces  a right 

A[z]-module  structure  on  A"  with  addition  and  scalar  multiplication 
given  by 


(v  + w)(k)  = v(k)  + w(k) 

v(  I z^a^.)  = I a^-Sp'(v) 

1 i 


(2.7) 
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where  v,w  e A^  a^-  e A and 

(SfvKk)  = F(k)(sJ-lv)(k-l)  for  1 = 12 
(Spv)(k)  = v(k). 

Let  p(2)  be  an  element  of  A[f]  given  by  Then  the  degree 

of  p(z),  denoted  by  deg  p(z),  is  the  maximum  va?ue  of  1 such  that 

ei  f 0.  If  deg  p(z)  = „ and  a„  = 1,  then  p(z)  is  said  to  be  monic.  Let 
A((rl))  denote  the  set  of  all  formal  Laurent  series  of  the  form 

00 

iLn  ^ " 0 Z . 

With  the  usual  addition,  and  with  multiplication  defined  by  (2.6)  and 

az->  = z-l(<,-la)  , a c A , 

A((z-1»  is  a noncommutative  ring  with  identity,  called  the  skew  ring  of 

formal  Laurent  series  over  A.  The  skew  ring  of  formal  power  series 
A[[z“^]]  defined  by 

A[[z-i]]  = {a  in  A((z-1))  : a = Y z‘'a  } 

i=o 

IS  a subnng  of  A((z-1)).  a polynomial  p(z)  e A[z]  is  said  to  be 

invertible  if  and  only  if  there  exists  a q(z)  in  A((z-1))  such  that 
P(z)q(z)  = q(z)p(z)  = 1. 
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A single-input  n-dimensi onal  system  (F,g)  is  said  to  be  N-cyclic 
with  N > n if  and  only  if  for  any  v e there  exist  ag,  aj,  . . . , 

afj_l  e A such  that 


N-1 

V = I a.Sp^g  . 

i=0  ^ 

From  the  multiplication  of  the  A[z]-module  structure  on  A^  defined  by 

(2.7),  we  can  express  the  above  equation  as  v = gir(z),  where  tt(z)  is  the 
polynomial  in  A[z]  given  by 

N-1 

tt(z)  = I z^a. 
i=0  ^ 

We  note  that  the  notion  of  N-cyclicity  is  equivalent  to  N-step 
reachability. 

Given  g e A*^,  let  Ann  g denote  the  annihilator  of  g defined  by 

Ann  g - {p(z)  e A[z]  : gp(z)  = 0).  Then  Ann  g is  a right  ideal  of  the 

ring  A[zL  KAMEN  and  HAFEZ  [1979]  have  shown  that  (F.q)  is  n-cvr.lir 

jn-step  reachable)  if  and  only  if  Ann  q ° afzlArzI.  where  a(z)  is  a 

monic  polynomial  of  degree  n.  This  monic  polynomial  a(z)  is  called  the 

9-  Let  a(z)  = z"  + z’a  From  (2.7),  we  see  that 

1=0 

ga(z)  = Sp'^g  + a.Sp''’g  . 

i=0  ^ ^ 


Therefore,  [ag  aj  . . . ~ -0^.^  ^(Sp^g),  where  the  prime  denotes 

transposition  (see  KAMEN  and  HAFEZ  [1979]). 
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From  the  results  of  KAMEN  and  HAFEZ  [1979],  if  a single-input 
system  (F,g)  is  n-step  reachable,  then  there  exists  a coordinate 
transformation  x(k)  = T(k)x(k)  such  that 

F(k)  = T(k+l)F(k)T-l(k)  , 

(2.8) 

g(k)  = T(k+l)g(k) 


where 


F(k)  = 

0 1 ...  0 

0 

• 

, g(i<)  = 

1 

• O 

I ...  * 1 

• 

• 

-ao(k)  -aj(k)  . . . -an_2(k) 

• 

. 1 . 

and  [ag  aj  . . . a^.i]'  = -C^-l(Sp''g)  . 

If  we  let  E,(k)  = [g(k)  (S.g)(k)  . . . (S.^-lg)(k)],  then  T(k)  is  given 
by  ^ 

(2.9)  T(k)  = CJk-l)Cr,-^(k-l)  . 

The  form  (2.8)  is  called  the  control  canonical  form  whirh  is  the 
discrete-time  counterpart  to  SILVERMAN  [1966]'s  phase  variable  form  in 
the  continuous-time  case. 


CHAPTER  THREE 

AUGMENTATION  OF  THE  SYSTEM  STATE 


In  this  chapter,  we  shall  present  a constructive  method  for 

augmenting  a N-step  reachable  time-varying  system  so  that  the  augmented 
system  has  a control  canonical  form. 

If  a given  n-dimensional  discrete-time  time-varying  system  is 
N-step  reachable  with  N > n,  we  will  first  show  that  it  is  possible  to 
turn  the  given  system  into  a system  having  a fixed  basis  bv  simply 
augmenting  the  system  state.  This  fixed  basis  leads  to  the  construction 
of  a control  canonical  form  for  the  augmented  system.  Then,  via  the 
control  canonical  form,  we  can  obtain  a feedback  gain  matrix  which 
results  in  the  closed-loop  system  matrix  of  the  augmented  system  being 
algebraically  equivalent  to  a constant  matrix  with  arbitrary  assignable 
eigenvalues.  Therefore,  by  using  dynamic  state  feedback,  which  consists 
of  the  augmentation  plus  the  feedback  gain  matrix,  we  have  an 
assignability  type  result  for  N-step  reachable  time-varying  systems. 

This  approach  to  assignability  yields  a controller  having  a much  lower 
order  in  general  than  the  controller  derived  by  KAMEN  et  al . [1985]. 

More  details  about  the  construction  of  feedback  control  systems  will  be 
given  in  Chapter  Four. 

Let  (F,G)  denote  the  m-input  n-dimensional  time-varying  system 
given  by  the  state  equation 


19 


20 


(3.1)  x(k+l)  = F(k)x(k)  + G(k)u(k)  . 

Given  a positive  integer  r,  let  D(k),  B(k),  E(k)  be  r x n,  r x r, 
time-varying  matrices,  and  consider  the  system 

(3.2)  w(k+l)  = D(k)x(k)  + B(k)w(k)  + E(k)u(k) 

The  system  given  by  (3.2)  is  an  augmentation  of  the  given  system  ( 
The  r X 1 vector  w(k)  is  the  state  vector  of  the  augmentation. 

Combining  the  system  equation  (3.1)  and  the  augmentation  (3.2 
obtain  the  following  state  equation  for  the  augmented  system 


(3.3) 


'x(k+l)' 

r — 

O 

IL. 

1 

'x(k)' 

’G(k)' 

= 

+ 

w(k+l ) 

1 

CO 

o 

1 

w(k) 

E(k) 

u(k) 


Figure  3.1  shows  the  block  diagram  of  the  augmented  system. 


/\  A 

In  the  following,  we  let  F(k),  G(k)  denote  the  matrices  defined  by 


F(k) 


’F(k) 

0 

. G(k)  = 

G(k)' 

,D(k) 

8(k). 

,E(k)_ 

r X m 


.1). 


, we 


To  derive  the  main  results,  we  need  the  following  proposition. 
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Given  System 


Fig.  3.1  The  block  diagram  of  the  augmented  system 
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(3.4)  PROPOSITION.  Let  X(L)  be  ae  „ x t matrix  over  A wUh  t > „. 

Suppose  that  ra„L  X(k)  = „ a„  integers  L.  Then  there  always  exists 
an  (A-n)  x X,  matrix  Y(k)  over  A such  that 


rank 


^ for  all  integers  k. 


proof.  Since  rank  X(k)  = n for  all  integers  k.  by 
version  of  Dolezal's  theorem  , there  exist  an  l x l 
^(k)  and  an  n x n invertible  matrix  W(k)  such  that 


the  discrete-time 
invertible  matrix 


X(k)M(k)  = [w(k)  [ 0] 

Let  V(k)  = [0  «i(l<)]M-l(k),  where  UjCk)  is  an  arbitrary 
X (^.-n)  matrix.  Then  we  have 


invertible 


(^-n) 


■ X(k)* 

= 

’ W(k) 

0 

. Y(k) 

• 

0 

Wi(k) 

M"l(k) 


Si  nee 


W(k) 

0 


is  invertible 


Wi(k) 


X(k)  is  also  invertible.  Therefore, 
Y(K) 


there  exists  a Y(k)  such  that 


is  not  unique. 


X(k) 

Y(K) 


is  invertible.  We  note  that  Y(k) 


□ 
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Now  we  have  the  following  results  in  the  single-input 


case. 


(3.5)  THOTM.  Suppose  that  (F.g)  is  N-step  reachable,  where  N > n. 
Then  there  always  exist  an  (N-n)  x n matrix  D(k).  an  (N-n)  x (N-n) 
matrix  B(k)  and  an  (N-n)  x 1 vector  e(k)  such  that  the  augmented 
N-dimensional  system  (F,g)  is  N-step  reachable,  where 


F = 


— 

F 

0 ' 

9 

g = 

— ■ 

g 

D 

B 

e 

P520F-  Since  (F,g)  is  N-step  reachable,  rank  CN(k)  = n for  all  integers 
k,  where  C(j(k)  = [g  Spg  . , . Sp"*g],  From  Proposition  3.4,  it  is 
always  possible  to  choose  (N-n)-element  column  vectors  eo(k),  ei(k), 

• • • . eM_2(k)  such  that 


rank 


9(K)  (Spg)(k)  . . , (ScN-lg)(k) 


ei(k) 


(k) 


for  all  integers  k, 


Let  C(k)  denote  the  N x N mat 


nx 


(3.6)  C(k)  = 


g(k)  (Spg)(k)  . . . (Sp^“lg)(k) 


ei(k)  . . . 
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Define 


e(k)  = eo(k)  . [D(k)  B(k)]  = [ei(k)  . . . e^,(k)](aC)-l  (k) , 

where  e^Ck)  Is  an  arbitrary  (N-n)-element  column  vector,  and  let 


F(k)  = 

'F(k) 

0 ' 

. g(k)  = 

■g(k)' 

_D(k) 

B(k)_ 

_e(k) 

By  applying  the  operator  S.  to  both  sides  of  (3.6),  we  have 


(S.C)(k) 

F 


(SF9)(k) 

ei(k) 


(SF^’^g)(k)  (SpNg)(k) 
eN-i(k)  e|^(k) 


Therefore,  C(k)  is  the  reachability  matrix  of  the  augmented  system 

A 

(F,g),  and  thus  the  theorem  is  proved.  [- 

Since  the  N-dimensional  augmented  system  is  N-step  reachable  , we 
can  apply  the  result  of  KAMEN  and  HAFEZ  [1979]  to  construct  a control 

canonical  form  for  the  augmented  system.  More  precisely,  we  have  the 
following  result. 


(3.7)  THEOREM.  If  (F,g)  is  N-step  reachable,  then  there  is  an 
N-dimens1onal  augmented  system  (F,g)  and  an  N x N invertible  matrix 
T(k)  such  that  T(ktl)F(k)T-'(k)  and  T(Ktl)g(k)  are  in  control 
canonical  form. 
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PROOF.  Combining  (2.8),  (2.9)  and  Theorem  3.5,  the  theorem  follows.  [] 

(3.8)  REMARK.  The  control  canonical  form  (F,g),  where  F = 

• A 

T(k+l)F(k)T(k)  and  g = T(k+1)  g(k),  is  called  a generalized  control 
canonical  form  in  the  single-input  time-varying  case.  [] 

We  illustrate  the  above  results  by  the  following  example. 

(3.9)  EXAMPLE.  Consider  the  single-input  time-varying  system  (F,g) 
given  by 


» 

1— » 

o 

1 

. g(i<)  = 

I 

1 

1 

_od  ( k ) 1 

1 

o 

1 

where  od(k)  = (1  , k is  odd 
(O  , k is  even 

As  in  Example  2.2,  the  reachability  matrix  is 


C3(k)  =11  1 


_0  od(k)  1 


and  this  system  is  3-step  reachable. 
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Let  C0(k)  = 


1 1 1 

0 od(k)  1 

1 0 1 


Then,  rank  C2(k)  = 3 for  all  integers  k 


Define  e(k)  = 1,  e3  (k)  = 0.  [D(k)  B(k)]  = [0  1 OlfaC’^k))  = [1  0 -1], 


Then  we  have 

tl 

t Ll_ 

1 

0 0 ' 

, g(k) 

= 

■ 1 ■ 

od(k)  1 0 

0 

1 

0 -1  , 

1 . 

Since  C"^(S.^g)  = 

F 

[1  -1 

-1]' 

, the  control  canonical 

form 

augmented  system 

is 

F(k)  = 

■ 0 1 

0 ' 

> 

g(k)  = 

■ 0 

0 0 

1 

0 

--1  1 

1 . 

1 

and  the  coordinate  transformation  is  T(k+1)  = 


where  ev(k)  ” ( ^ s k is  even 
V 0 , k is  odd. 


"O 

0 

1 

1 

1 

od(k) 

ev(k) 

1 

-ev(k) 

ev(k) 

1 

od(k) 

[] 


Now  let  us  consider  the  m-input  n-dimensional  time-varying  system 
(F,G).  In  the  following,  we  assume  that  rank  G(k)  = m for  at  least  one 
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integer  value  of  k.  If  rank  G(k)  < m for  all  integers  k,  then  the 
augmentation  will  have  the  structure  shown  in  Fig.  3.2,  where  u is  the 
input  of  the  system  (F,G)  and  G,  Mi(k)  are  defined  in  Chapter  Two. 

(3.10)  HEOREM.  Suppose  that  (F,G)  is  N-step  reachable  and  rank  G(k)=m 
for  at  least  one  integer  value  of  k.  Then  for  some  integer  q where  N < 
q < mN  there  exists  a (q-n)-dimensional  augmentation 

(3.11)  w(k+l)  = D(k)x(k)  + R(k)w(k)  + E(k)u(k) 

such  that  the  q-dimensional  augmented  system  (^G)  is  q-step  reachable 
with  a fixed  basis. 


proof.  We  would  like  to  prove  the  above  theorem  step  by  step  in  a 
constructive  way. 


STEP  I.  Using  the  selection  process  described  in  Chapter  Two,  determine 

the  generalized  controllability  indices  of  the  system,  nj,  

n^.  For  1 < i < m , we  can  then  express  Sp^^’g^  as  an  A-1  inear 
combination  of  previously  selected  vectors;  that  is 


(3.12) 


Sp^g 


= - I 

j=0 


m 

I 

l=l 


r^VsJ 

j ^ F 


i=l,  2, 
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Input 

Transformation  Given  System 


Fig.  3.2  Augmented  system  when  rank  G(k)  < m 
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where  e A with 

J 


= 0 
m 


for 


^ = 1,  i + 1, 


. . , m 


r^T  = = 

^ * 


r^^  = 0 


for  £ = 1,  2, 


, m 


STEP  II.  Choose  q vectors  e,j(k),  1 < j < m,  0 < 1 < nj  -1,  where 
eij(k)  are  (q-n)  x 1 vectors,  such  that 


rank 


9l  Spgi  . . . Sp"'  9i  92  Spgj  . . . g,„  ...  SF""''g„ 


L^Ol  ni  • . . eo2  ejj  . . . eq^  . . . 0(„  i)„ 


at  all  times. 


This  is  always  possible  from  Proposition  3.4.  Let  Cq(k)  denote  the 
q X q matrix 


(3.13)  Cq 


9l  Spgj  . . . Sp  * 9i  ...  ...  SF""'*g„- 

,®01  eil  ...  e(„j.i,i  eo2  . . . eo„  . . . 
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Using  the  coefficients  in  (3.12),  calculate  e„  . , i = 1 2 

j ri  ^ I * 9 9 

such  that 


• 5 TTl 


1 m 


.n'i 


® " > i ~ 1 , 2 , , 


, m 


(3-14)  e i = - I ^ 

^ j=0  1=1  J 

STEP  III.  Define  E(k)  = [eoi(k)  eQ2(k)  . . . eom(k)]  and 


[D(k)  B(k)]  = X(k)(aCq)“^(k) 


where  X(k)  = [e„  . . . ejj  . . . . . . e„ 


Let  F = 


F 

0 

and 

A 

G = 

’ g' 

_D 

B _ 

E 

Apply  the  operator  to  both 

F 


sides  of  (3.13).  Then  we  have 


S.C, 


SpOi  . . . Sp  Sp92  . . . Spg^  . . . Sp\^ 


®11  • • • ei2 

From  (3.12)  and  (3.14),  we  see  that 


®lrn  • • • 


n. 

1 m 


^F  ’ jL  Jl  ’ ■*  = 2,  . . .,  m 
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where  g.  is  the  column  vector  of  G and  g S.^g  are  given  by 

F ^ 


9i  = 

'si  ' 

F 

. . 

Therefore,  Cq  will  be  a part  of  the  reachability  matrix  for  the 
augmented  system  (F,G),  where  Cj^  is  given  by 

c,^  = [G  si  . . . Si"i]  . 

F F 

Then  the  theorem  follows.  j--j 

We  should  note  that  by  equation  (3.14),  the  controllability  indices 
of  the  augmented  system  are  the  same  as  the  indices  of  the  given  system 
which  are  determined  in  Step  I of  the  proof  of  Theorem  3.10. 

We  now  define  a control  canonical  form  for  the  time-varying 
system.  Suppose  that  the  m-input  n-dimensional  time-varying  system 
(F,G)  has  a coordinate  transformation  x(k)  = T(k)x(k)  such  that 
F(k)  . T(kH)F(k)T-l(k)  and  G(k)  = T(k+l)G(k)  are  in  the  following  form 


■^11 

''12  • • . Fj^ 

» 

G = 

■ ■ 

''21 

• * • ''2m 

G2 

• 

« 

• 

• 

• 

• 

iml 

• * • ''mm  _ 

- '^m  . 
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where 

^ii  = 

°(rii-l)xl  ^n^--l 

ifj 

*^{n^--l  )xnj 

[x  1 X ...  X 

X ...  X 

Gi  = 

°n^-  x(i-l) 

°(n.j-l)x(n-i+l) 

1 X ...  X 

and  X's  are  coefficients  in  A. 


Then  the  form  (3.15)  is  called  the  control  canonical  form  nf  the 
multivariable  time-varying  system.  For  the  control  canonical  form  of 
the  time-invariant  case,  the  reader  may  refer  to  POPOV  [1972]  .KAILATH 
[1980],  and  CHAN  and  WANG  [1978]. 

Since  the  augmented  system  constructed  in  the  proof  of  Theorem  3.10 

has  a fixed  basis,  we  can  construct  the  control  canonical  form  for  the 
augmented  system. 

(3.16)  theorem.  Let  (F,G)  be  the  augmented  system  constructed  in  the 
proof  of  Theorem  3.10.  Then  there  is  a coordinate  transformation  x(k)  = 
T(k)x(k)  such  that  F(k)  = T(ktl)f(k)T-l(k)  and  G(k)  = T(kH)G(k)  are  in 
control  canonical  form. 

proof.  This  proof  is  a generalization  of  the  proof  in  the  time- 
invariant  case.  The  reader  may  refer  to  CHAN  and  WANG  [1978]. 
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us  denote  G [g^ , , g^].  Using  the  coefficients 

equation  (3.12),  define  the  following  q-element  column  vectors 


in 


F j=0  £=1 


n,--J  ' r 


i 1.2,  . . . m ; t = 0,1,  . . . ,n^ 


Let  be  the  q x q matrix  defined  by 


vm 

""m-l 


Then  n can  be  shown  that  v'T,  = 0 if  and  only  If  , = 0.  Hence,  T,  is 
invertible  (see  POPOV  [1972]).  Now.  we  must  show  that  T,  is  a 
coordinate  transformation  such  that 

(3.19)  T F = FaT 

^ X 

(3.20)  T^G  = G 

where  F and  G are  in  control  canonical  form.  From  (3.17)  and  the 
coefficients  in  equation  (3.12),  we  have 


><0  = 9i  ^ 

t=l  ' 


1.  2,  . . ., 


m 


If  we  write  these  equations  in  a matrix  format,  we  have 
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(3.21) 

lxq  . . 

1 — 1 
E O 
X 

• 

= Cgj  . . . 

1 

0 

■^2  12 
^2 

1 

• • 

■n  1 

m Im 
. a r 

n 

m 

• 

• 

• 

• 

• 

• 

• 

• 

^■'’m^(ni-l)in 

1 

0 

• 

• • 

Let 

Q = 

1 

0 

-ng  12  -n 

^ r'  • • • (7 

"2 

1 

m Im ' 
tn 

• 

• 

• 

-1 

•l  b'2 
0 1 

• • • u 

k2m 

• • • IJ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

l^(m-l)m 

. 0 

1 

1 

Then  we  have  T^^G  = G by  inserting  zero  rows  in  Q,  where 


(3.22)  G 
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To  prove  (3.19)  , we  apply  S.  to  both  sides  of  (3.17).  Then 


(3.23)  S xj  = S^‘g  t J f 

^ F j=0  £=1  n.-j  p 


where 


1 ? , -{n,-t-D 

- Vl  ■ J.  (“ 


£=1 


i r"  ",  -1  -(n.-t-l) 

= Vl  - f9i  • • . g^](Q  ^)(Q)  ' 


^t+1  LXq  . . . XnJ 


, li 
n^.-t-l 


r,li 


-(n.-t-l)  . 

„ 1 mi 

a r 

n^. -t-1 


n.j-(t+l) 

• 

: 

L‘'n.-(t+l)  J 


i = 1,  2, 


n.j-t-1 


. 1 

n^.  -t-1' 


. , m ; t - 0,  1,  . . , j n.j-1. 


(Q) 


-(r,i-t-l)  j, 


n^.-t-l 


•(n.-t-l) 

1 ^ mi 


n^-t-1 


By  inserting  zeros  in  (3.23)  and  writing  these  equations  in  a matrix 
format,  we  have 


= T F 

n X X 
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where 


(3.24)  F. 


1 1 


-a 


n 


= 


■aJJ  -a}J 


n.-l 


^(n^--l)x(m-ixl) 

• • • 

T » J > = 1 , 2 , . . . , m , i j . 


Letting  T(k)  - we  have  T(k)F(k)  = T(k+l)F(k)  and  G(k) 

T(k+l)G(k). 


□ 


The  control  canonical  form  (F.G)  for  the  augmented  system  (F.G)  is 

A generalized  control  canonical  form  in  the  rwlti variable  time- 

varying  case.  We  now  illustrate  the  above  results  by  the  following 
exampl e. 


37 


(3.25)  EXAMPLE.  Consider  the  system  (F,G)  given  by 


F(k)  = 

’ 1 

0 ' 

, G(k)  = 

’ 1 1 ' 

O 

1 

1 

_ 0 od(k) 

where  od(k)  = H , k is  odd 

( 0 , k is  even 
Then  C2(k)  is  given  by 


ev(k)  = ( 1 

( 0 


k is  even 
k is  odd 


C2(k) 


1111 
. 0 od(k)  0 ev(k) 


which  has  rank  2 for  all  Integers  k.  Thus  the  system  is  reachable  in  2 

steps.  But  C2(k)  does  not  have  a fixed  basis.  Now,  let  us  construct 
the  augmented  system. 


iIEP_J_.  We  select  g^,  g2  and  Spg2  from  C2(k)  and  express  Sp^g2  and  Spg^ 
as  A-linear  combinations  of  these  previously  selected  vectors.  Then 

(SFgi)(k)  = gi(k) 

(Sp^92)(k)  = g2(k) 

The  generalized  controllability  indices  are  n^  = 1 , n2  = 2 and  the 
dimension  of  the  augmented  system  is  q = 3. 
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SJEPJl,  If  we  choose  eQi  = 0,  eQ2  = -1,  e^  = 1, 


Cs(k) 


1 1 1 

0 od(k)  ev(k) 
0 -1  1 


and  rank  C3  3 at  all  times.  Now,  determine  e^^  and  e22  from  Step  I 
and  egi,  eg2,  e^j.  Then  we  have 

eii  = 0 

622  = -1 

STEP  HI.  Define  E(k)  = [0  -1]  and 


[ D B ] = X(k)(aC(k))"l 
= [0  1 -1  ] 

■ 1 

-2 

od(k)-ev(k) 

0 

1 

-od(k) 

0 

1 

ev(k ) 

=[00-1]. 

Hence,  the  augmented  system  (F,G)  will  be 


*1 

0 

1 — 

0 

A 

. G = 

‘1  1 

0 

1 

0 

0 od(k) 

_0 

0 

-1  _ 

\ 

1 

0 

1 
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Now,  let  us  construct  a generalized  control  canonical  form.  From  Step 
II. 

(S^9l)  = Qi  , 

F ^ ^ 

(S?92)  = 92 

»0  = 9l  . xi  = S.gj  , = 92  , 

and  let  = [xg  Xg]  . Then  we  have 

= [xj  xj  xf]  = T/ 

G = T,G 

Therefore,  the  generalized  control  canonical  form  is 


II 

1 Li. 

1 

0 

o 

J 

G = 

‘ 1 

0 ' 

0 

0 

1 

0 

0 

_ 0 

1 

0. 

0 

1 

In  the  following  chapter,  we  shall  construct  a feedback  gain  matrix 
for  the  generalized  control  canonical  form  which  results  in  the  closed- 
loop  system  matrix  being  algebraically  equivalent  to  a constant  matrix 
with  arbitrary  assignable  eigenvalues.  Later,  the  augmentation  approach 
will  be  applied  to  the  design  of  output  feedback  controllers. 


CHAPTER  FOUR 
FEEDBACK  CONTROL 


This  chapter  is  divided  into  two  parts.  In  the  first  part,  we 
shall  show  how  a feedback  gain  matrix  can  be  constructed  for  augmented 
N-step  reachable  time-varying  systems  using  the  generalized  control 
canonical  form,  resulting  in  assignability  of  the  closed-loop  system 
dynamics.  If  an  n-dimensional  time-varying  system  is  N-step  reachable 
with  N > n,  then  it  is  possible  to  assign  the  coefficients  of  the 
closed-loop  characteristic  polynomial  by  employing  a dynamic  controller 
consisting  of  the  feedback  gain  matrix  and  the  augmentation.  The  system 
matrix  of  the  resulting  closed-loop  system  is  algebraically  equivalent 
to  a constant  matrix  with  arbitrary  assignable  eigenvalues.  This 
assignability  result  for  N-step  reachable  time-varying  systems  was  first 
obtained  by  POOLLA  [1984]  and  KAMEN  et  al.  [1985].  But,  we  believe  that 
our  scheme  for  finding  a controller  is  much  more  flexible,  and  our 

controller  has  a much  lower  order  in  general  than  the  controller  derived 
by  KAMEN  et  al . [1985]. 

In  the  second  part  of  this  chapter,  the  augmentation  approach  is 
applied  to  design  output  feedback  controllers.  In  this  case,  we  can 
also  get  the  coefficient  assignability  result.  Later,  the  separation 
property  of  the  observer-controller  design  procedure  will  be  proved. 
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The  augmented  system  constructed  in  Chapter  Three  for  the  N-step 
reachable  time-varying  system  (F,G)  is  given  by  the  following  state 
equation 

(4.1)  x(k+l)  = F(k)x(k)  + G(k)u(k) 

where  x(k)  - x(k)  , F(k),  G(k)  are  q x q , q x m matrices  such  that 

w(k) 


1 

n 

o 

J 

A 

, G = 

1 

CD 

L. 

1 

cc 

Q 

1 

1 

LlJ 

1 

The  problem  is  to  find  a nonsingular  m x m matrix  P(k)  and  an  m x q gain 

A 

matrix  L(k)  such  that 

(4-2)  u(k)  = P(k)v(k)  - L(k)x(k) 

makes  the  closed-loop  state  equation 

x(k+l)  = (F(k)  - G(k)L(k))x(k)  + G(k)P(k)v(k) 

have  the  preassigned  characteristic  polynomial  where  v(k)  is  a reference 
input. 

From  Theorem  3.16,  there  exists  a coordinate  transformation  x(k)  = 

A 

T(k)x(k)  such  that  under  this  transformation,  the  state  equation  (4.3) 
becomes 


(4.4) 


x(k+l)  - (F(k)  - G(k)L(k))x(k)  + G(k)P(k)v(k) 
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where 

F(k)  = T(k+l)F(k)T~^(k)  , G(k)  = T(k+l)G(k) 

and  L(k)  = L(k)T"l(k) 


The  matrices  F,G  have  the  generalized  control  canonical  form 
(3.15). 

Let  us  determine  the  matrix  P(k)  such  that 


(4.5)  GP  = block  diag  {[0 


0 1]' , n^-  X 1 , i=l. 


m } 


where  the  above  notation  denotes  the  matrix  whose  diagonal  elements  are 
riT  X 1 column  vector  [0  . . . 0 1]'  for  1 < i < m.  Then  from  (3.22), 
P(k)  has  the  form 


P 


1 

0 


p(m-l)m 

1 


where  are  the  values  in  A determined  from  the  elements  of  the  matrix 
G.  Define 
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where 


^•i  = 

L-ag 

- fg 

• "^n]-l  ■ 

1 — 1 

1 

•r-  C 

= [■ 

■4' 

• • • 

-4\.o 

• 

J 

Here,  fj^' 

, 1 < i 

< m. 

0 < j 

< are 

the  preass 

igned  coefficients 

of 

the  system  matri 

X of 

equation  (4.4)  and 

ajJ,  1 < i 

,j  < m,  0 < £ < 

1 , 

are  the  elements 

of  the  matrix 

F.  Then 

we  have 

(4.6)  F - 

GPL  = 

block 

diag 

{ 

'O  1 0 

• 

...  0 
• 

, n^-xn^-,i=l, 

• • • » 

m} 

• 

• 

• 

1 

^0  . . 

-fii 

Therefore 

, if  we 

let 

'(k)  = 

P 

(k)L(k). 

the  desired 

feedback  gain 

mat 

rix 

is 

(4.7) 

L(k) 

= P(k)L(k)T(k) 

Figure  4.1  shows  the  block  diagram  of  the  feedback  system  structure 
for  the  augmented  system.  The  dynamic  controller  for  the  given  system 
(F,G)  is  the  augmentation  plus  the  feedback  gain  matrix  L(k).  This 
controller  will  results  in  the  augmented  system  having  a preassigned 
closed-loop  system  dynamic. 
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Fig.  4.1  Feedback  system  structu 
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Wg  illustratG  the  abovG  rGSult  by  the  following  GxamplG. 


(4.8)  EXAMPLE.  Let  us  consider  the  augmented  system  constructed  in 
Example  3.25.  The  augmented  system  (F,G)  is  given  by 


' 1 

0 

0 ‘ 

. G(k)  = 

■ 1 

1 ■ 

0 

1 

0 

0 

od(k) 

_ 0 

0 

-1  _ 

_ 0 

-1 

and  the  control  canonical  form  (F,G)  for  the  augmented  system  is 


F = 

’ 1 0 0 ■ 

> G = 

* 1 o“ 

0 0 1 

0 0 

1 

o 

1—* 

o 

1 

o 

1 

where  the  coordinate  transformation  T(k)  is 


T(k) 


1 

0 

0 


-2  od(k)  - ev(k) 

1 ev(k) 

1 - od(k) 


Suppose  that  we  want  the  closed-loop  eigenvalues  at  1/2,  0,  1/4.  Then 


L 


1/2  0 0 

0 1 1/4 


In  this  case,  P = I3  and 
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L(k) 


1/2 

0 


-1  l/2{od(k)  - ev(k)>  ■ 

5/4  ev(k)  - 1/4  od(k) 


Figure  4.2  shows 
0 and  the  initial 


the  simulation  results  of  the  above  example  when  v(k) 


conditions  are  x(0)  = 


2 

1 


. w(0)  =0.2  . 


□ 


We  now  consider  the  MIMO  linear  time-varying  discrete-time  system 
given  by  the  input/output  difference  equation 

Pi"!  P2-I 

(4.9)  y(k+pi)  + I Qi(k)y(k+i)  = I R.-(k)u(k+i) 

i=0  i=o 


where  pj  > P2,  are  p x p,  and  are  p x m time-varying  matrices. 
Equation  (4.9)  can  be  rewritten  in  the  following  form; 


Pl-1 

(4.10)  y(k+l)  = I Q.(k-pj+l)y(k-p,+l+i) 

i=0 

P2-I 

I Pi  (k-pi+l)u(k-pi+l+i ) 
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a)  Input  u(k)  = [u^(h)  U2(k)]' 


\ 

\ 

\ 

L 

\ 

*1 

\, 

: ') 

^ 

* » » ‘ * * * M 

> 1 1 M >, 

-Li  »>»«>»  » 

i 

* > * * ‘ ‘ ■ 

5 

^ * ^ » * ‘ * ‘ I 

2.«  3.0  4.0  5.0 


b)  States  x(k)  = [Xj(k)  X2(k)]'  and  w(k) 


Fig.  4.2  Simulation  results  of  Example  4.8  (v(k)  = 0.  x (0)  = 2 
X2(0)  = 2,  w(0)  = 0.2)  > J 
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The  above  equation  can  be  expressed  as  follows: 


y(k+l)  - -Qp^_i(k-p;^+l)y(k)  + r^Ck) 

r^fk+l)  = -Qp^_2(k-pi+2)y(k)  + r2(k) 

• • 

• • 

'•pj.pj.llM)  = -Qp2(k-P2)y(k)  + rpj.p^(^, 

' -0p2(l<-P2+l)y(k)  + rp^.i(k-p2+l)u(k) 

* '■pi-P2+l<^> 

• • 

• • 

• • 

^pj-l(l<+l)  = -QoO<)y(k)  + ro(k)u(k) 


If  we  define  the  state  vector  x'(k)  = [rj;^_i(k)  . . . r'(k)  y'(k)]', 
then  we  have  the  state  space  model 


x(k+l)  = F(k)x(k)  + G(k)u(k) 

\ 4 , 1 i j 

y(k)  = H(k)x(k) 

where 


4.12a) 


F(k) 


0 


0 -0o(k) 

0 -Qi(k-l) 
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(4.12b)  G(k) 


Ro(k) 

Rl(k-l) 

®pxm 

®pxm 


, H = [0, 


• • ^p  ^p^ 


Let  n = pjp.  Then  F is  n x n , G is  n x m , and  H is  p x m.  We  note 
that  the  above  realization  is  observable  in  n steps  at  all  times.  The 
observable  realization  (4.12a,b)  is  a discrete-time  time-varying  version 
of  Fuhrmann's  realization  (see  FUHRMANN  [1976]  and  POOLLA  [1984]). 

If  the  n-dimensional  observable  realization  (4.12  a,b)  is  N-step 
reachable  with  N > n , we  can  use  the  augmentation  approach  to  obtain  a 
combined  observer-control ler  compensator  which  assigns  the  coefficients 
of  the  characteri Stic  polynomial  of  the  closed-loop  system 
arbitrarily. 


(4.13)  THEOREM.  If  the  m-input  n-dimensional  time-varying  system 
(F,G,H)  defined  by  (4.12a,b)  is  N-step  reachable  for  some  integer  N > n 
there  is  an  output  feedback  compensator  of  dimension  q,  where 
N < q < mN,  that  results  in  the  closed-loop  system  matrix  being 
algebraically  equivalent  to  a constant  matrix  with  q arbitrary 
assignable  eigenvalues. 
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proof.  Given  the  realization  (4.12a,b),  define 


M(k) 


-Qo(k) 

-Ql(k) 


Qpj-l(k-Pi  +1) 


and  consider  the  observer 

(4.14)  Xo(k+l)  = F(k)xo(k)  + M(k)(y(k)  - Hxo(k))  + G(k)u(k) 

If  (F,G)  is  N-step  reachable,  for  some  integer  q where  N < q < mN  there 
exists  a (q-n)-dimensional  augmentation 

w(k+l)  = D(k)xo(k)  + 8(k)w(k)  + E(k)u(k) 

where  D(k),  B(k),  E(k)  are  determined  by  the  method  described  in  Chapter 

A A 

Three  so  that  (F(k),G(k))  has  a fixed  basis.  From  Theorem  3.16,  there 
exists  a coordinate  transformation  T(k)  such  that  T(k+l)F(k)T'^(k) 

A 

T(k+l)G(k)  are  in  control  canonical  form.  Then  we  can  choose  an  m x q 
feedback  matrix  L(k)  so  that  T(k+l)(F(k)  - G(k)L(k))T-l(k)  has  the 
preassigned  matrix  of  the  form  (4.6). 
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Finally,  consider  the  feedback  control  law 


(4.16) 


u(k) 


-L(k) 


xo(k) 

w(k) 


Let  L(k)  CL2(k)  L2(k)],  where  Lj(k)  is  an  m x n matrix  and  L2(k)  is 
an  m X (q-n)  matrix.  Combining  (4.9),  (4.14),  (4.15)  and  (4.16),  we 
have  the  closed-loop  system 


where 


^CL(*^)''(k) 

xcL(k)  = 


■ x(k) 

■ G(k)' 

w(k) 

E(k) 

1 

O 

X 

1 

. G(k). 

Fa(k)  = 


F 

0 

-F+MH 


-GL^ 


B-ELc 


-GLi 

D-ELi 

F-MH 


Define  xcl(I<)  = P"^X0i_(k),  where 


)-l  = 


I 

0 

-I 


n 


n 


q-n 
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Then  we  have 


where 


Xcl(I<+1)  - FQL(l<)Xf;L(k)  + GQ|_(k)v(k) 


= P-^clP  - 


F-GL 


1 


D-EL 

0 


1 


-GLc 


B-ELc 


-GLi 


D-ELi 


F-MH 


G(k) 

E(k) 


2G(k) 


By  the  discrete-time  version  of  the  Lyapunov  transformation,  is 
algebraically  equivalent  to  x^l.  Since  F-MH  is  a fixed  matrix  and 


A A A 

F - GL 


F-GLi 

D-ELi 


-GL2 

B-EL2 


can  have  q arbitrary  eigenvalues,  the  closed-loop  system  matrix  FQ|_(k) 
is  equivalent  to  a time-invariant  matrix  with  q assignable  eigenvalues. 

□ 
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We  now  consider  the  SISO  case.  The  input/output  difference 
equation  is  given  by 


n-1  n-1 

(4.17)  y(k+n)  + I ai(k)y(k+i)  = I b.(k)u(k+i)  . 

i=0  i=o 


By  a similar  derivation  as  the  MIMO  case,  we  can  obtain  the  observable 
realization  (F,g,h)  of  (4.17)  such  that 


(4.18a)  F(k) 


. . -ao(k) 

. 0 -ai(k-l) 


^ ■3n-i(k-n+l) 


(4.18b) 


g(k) 


bo(k) 

bi(k-l) 

bn_i(k-n+l) 


. h = [0  ...  0 1] 


Then  we  have  the  following  theorem  which  is  analogous  to  Theorem  4.13, 


(4.19)  theorem.  If  the  system  (F,g,h)  defined  by  (4.18a,b)  is  N-step 
reachable  for  some  integer  N > n , there  is  an  output  feedback 
compensator  of  dimension  N that  results  in  the  closed-loop  system  matrix 
being  algebraically  equivalent  to  a constant  matrix  with  N arbitrary 
assignable  eigenvalues.  n 
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We  have  so  far  considered  the  construction  of  output  feedback 
compensators  for  the  time-varying  systems  given  by  the  input/output 
difference  equations  (4.9)  and  (4.17).  Now  let  us  consider  any  time- 
varying  system  that  is  observable  and  reachable  in  finite  steps.  In  the 
following,  we  prove  the  separation  property  of  the  observer-controller 
design  procedure  when  the  augmentation  approach  is  used. 

Let  the  system  (F,G,H)  be  given  by  the  state  equation 

x(k+l)  = F(k)x(k)  + G(k)u(k) 

(4.20) 

y(k)  = H(k)x(k) 

where  x e y e AP,  u e A'^.  Suppose  that  the  system  (F,G,H)  is 
Nj^-step  reachable  and  N2-step  observable,  where  Ni,N2  > n.  Then  by  the 
previous  results  there  exists  an  augmentation  for  (F,G)  such  that  for 
some  integer  q^  with  < q^  < mN^ , the  qj-dimensional  augmented  system 

^ A 

(Fj,  G)  is  q-step  reachable  and  has  a fixed  basis,  where 


= 

' F 

1 

o 

» 

A 

G 

G 

. 

1 

m 

1 

Let  L(k)  - [Lj(k)  L2(k)  ] be  the  feedback  gain  matrix  for  the  augmented 

A A 

system  (F,G).  Then  the  closed-loop  system  matrix  will  have  the  form 
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(4.21)  - GL 


F - GL^ 

. Di  - EiLi  B1-E1L2 


By  the  duality  principle,  there  also  exists  an  augmentation  for  (F,H) 
such  that  for  some  integer  P2  with  N2  < q2  < pN2,  the  q2-dimensional 

A A 

augmented  system  (F^.H)  is  q2-step  observable  and  has  a fixed  basis, 
where 


F 

0 


[H  E2] 


B2  J 


Consider  the  observer 


/V  A 


(4.22)  XQ(k+l)  = F2(k)xg(k)  + M(k)(y(k)  - H(k)x.(k))  + 


G(k) 

0 


u(k) 


where  XQ(k)  = 


XQ(k) 


is  the  state  of  the  observer,  and  M(k)  = 


Ml(k) 

M2(k)_ 


IS 


the  feedback  matrix  of  the  observer.  Then  the  closed-loop  system  matrix 
of  the  observer  will  be 


(4.23) 


MH 


F - MjH 
-MH2 


D2  - M1E2 
B2  - M2E2 


5 6 

From  the  previous  results,  we  have  the  augmentation  and  the 
feedback  control  law  such  that 


w(k+l)  = Di(k)xo(k)  + Bi(k)w(k)  + Ei(k)u(k) 

(4.24) 

u(k)  = -l-i(k)xo(k)  - L2(k)w(k)  + v(k) 

Combining  (4.20)  ,(4.22)  and  (4.24),  we  have  the  closed-loop  system 

(4.25)  XcL(k+l)  = FcL(k)xcL(k)  + GcL(k)v(k) 
where 


XQL(k+l)  = 

‘x(k+l) 

’G(k)  ■ 

w(k+l) 

Ei(k) 

xo(k+l) 

G(k) 

+ 

0 

1 X 
1 

. 0 . 

'^a(k)  = 


F 

0 


MjH 

M2H 


-GL2 

B1-E1L2 

-GL2 

0 


-GLi 

Di-EiLi 

F-MjH-GLi 

-M2H 


0 

0 


D2“M]^E2 

B2“M2E2 


Define  XQ|_(k)  = P"^xql(I^).  where 


(4.26) 


0 

■L 


I 


0 

qi-n 

0 

0 


0 

0 

L 


^2 


-n 
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Then  we  have 


(4.27) 


XcL(k+l)  = FcL(k)xcL(k)  + GcL(k)u(k) 


where 


Fcl  = 


Di  - E^Li 
0 
0 


-GLi 

B1-E1L2 

0 

0 


-GLi 

Dl 

F-M^H 

-M2H 


0 

0 

D2“M2E2 

B2~^2^2 


G = 


G 

El 

2G 

0 


By  the  discrete-time  version  of  the  Lyapunov  transformation , is 
algebraically  equivalent  to  Xql  and  from  (4. 21), (4. 23)  and  (4.27)  we  see 
that  the  controller  and  the  observer  can  be  designed  separately.  This 
is  the  separation  property  of  the  observer-control ler  design  procedure. 
Figure  4.3  shows  the  structure  of  the  combined  observer-control ler 
feedback  system. 
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Fig.  4.3  Combined  observer-controller  structure 


CHAPTER  FIVE 
CYCLIZABILITY 


In  this  chapter,  we  shall  use  the  augmentation  approach  to  derive  a 
different  constructive  method  of  the  design  of  feedback  control  systems 
for  N-step  reachable  time-varying  systems  from  Chapter  Four.  Since  the 
augmented  system  constructed  in  Chapter  Three  has  a fixed  basis,  it  was 
possible  to  obtain  a generalized  control  canonical  form  which  led  to  the 
design  of  a feedback  control  system  which  results  in  the  closed-loop 
system  matrix  being  algebraically  equivalent  to  a constant  matrix  with 
arbitrary  assignable  eigenvalues.  However,  it  may  be  quite  laborious  to 
transform  the  system  to  the  control  canonical  form.  It  is  therefore  of 
interest  to  obtain  a more  direct  method. 

If  a time-invariant  multivariable  system  is  reachable,  then  there 
exists  a feedback  gain  matrix  so  that  the  closed-loop  system  is 
reachable  through  one  of  the  inputs.  Therefore,  it  is  possible  to 
construct  feedback  control  systems  by  using  single-input  methods  without 
having  to  consider  the  multivariable  control  canonical  form.  Moreover, 
we  can  achieve  the  coefficient  assignability  by  this  method.  This 
notion,  which  is  called  cycl izabi 1 ity , was  first  introduced  by  HEYMANN 
[1968].  It  is  well  known  that  a time-invariant  system  is  reachable  if 
and  only  if  it  is  cyclizable  (see  HEYMANN  [1968]  and  HAUTUS  [1977]). 
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However,  in  the  time-varying  case,  the  notion  of  cyclizability  has 
not  been  thoroughly  investigated,  and  there  are  still  several  open 
questions.  For  example,  if  a time-varying  system  is  cyclizable,  then  it 
is  always  reachable.  But,  it  is  not  known  whether  or  not  reachability 
implies  cyclizability. 

To  study  the  time-varying  case,  we  begin  with  the  time-varying 
version  of  the  definition  of  cyclizability. 

An  m-input  n-dimensional  time-varying  system  (F,G)  is  said  to  be 
N-cyclizable  with  N > n if  there  exist  a matrix  Lqg  and  a vector 
ug  e A'^  such  that  (F-GLg,GuQ)  is  N-cyclic.  If  (F-GLg,  Gug)  is  n-cyclic, 
then  the  system  is  said  to  be  n-cyclizable.  We  note  that  there  is  no 
notion  of  N-cycl izabi 1 ity  in  the  time-invariant  case.  In  the  following 
lemma,  KAMEN  and  HAFEZ  [1979]  have  obtained  a necessary  and  sufficient 
condition  for  n-cycl izabi 1 ity  in  the  time-varying  case. 

(5.1)  LEMMA.  There  exist  a matrix  Lg  e A^^^  and  a vector  ug  e A>^  such 
that  (F-GLg,  Gug)  is  n-cyclic  (or  n-step  reachable)  if  and  only  if  there 
exist  Ug,  uj,  . . . , Up_;^  e A^  such  that 

rank  [Cg  • • • ?n-l^  ~ at  all  times 

where 

5g  = Gug  , 

= F(<?Ci_i)  + Gu^-  , i=l,2,  . . . , n-1.  [] 
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However,  in  the  time-varyi ng  case,  n-step  reachability  of  the 
system  (F,G)  does  not  always  imply  that  (F,G)  is  n-cyclizable  as  we  see 
in  the  following  lemma. 


(5.2)  LEMMA . Suppose  that  the  system  (F,G)  is  n-step  reachable.  If 
G(k)  - 0 for  at  least  one  value  of  k,  then  (F,G)  is  not  n-cyclizable. 

_PR00F_.  For  the  system  (F,G)  to  be  n-cyclizable,  there  must  exist  a 
matrix  Lq  e and  a vector  ug  e such  that  (F-GLg,  Gug)  is 
n-cyclic.  If  G(kj)  = 0 for  some  integer  k^ , then  G(ki)ug  = 0 for  any 
vector  ug  e A"^.  Therefore,  for  the  vectors  defined  in  Lemma  5.1, 

rank  [?g  ^ r\  at  time  kj,  and  thus  the  system  (F,G)  is  not 

n-cyclizable.  rg 


However,  even  if  (F,G)  is  not  n-cyclizable,  it  may  be  N-cyclizable 
for  some  finite  integer  N > n.  In  other  words,  there  may  exist  an  Lg 
and  an  ug  such  that  (F-GLg, Gug)  is  N-cyclic  when  the  system  is  n-step 
reachable. 


(5.3)  EXAMPLE.  Let  F(k)= 


’ 1 

0 ' 

and 

M 

CD 

k 

0 ' 

_ 0 

1 _ 

_ 0 

k_ 

Then  C G SpG  ] = 


k 0 k-1 

0 k 0 


0 

k-1 


which  has  rank  2 for  integers  k, 
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Thus  the  system  is  2-step  reachable.  But,  since  G(k)  = 0 at  k = 0,  by 
Lemma  5,2,  (F,G)  is  not  2-cycl izable.  However,  if  we  choose 


'-0  = 

1 

o 

o 

J 

and 

% 

= 

1 

. ^ 

_ 0 _ 

1 0 

and 

Gug 

* k ' 

• 

^-k  1 _ 

0 

The  system  (F-GLq,  Guq)  is  3-step  reachable.  Therefore,  the  system 
(F,G)  is  3-cycl izable.  [] 

Now  suppose  that  the  system  (F,G)  is  N-step  reachable  and 
M-cycl i zabl e,  where  M > N > n.  That  is,  let  us  suppose  that  there  exist 
an  Lg  and  an  ug  such  that  (F-GLg,  Gug)  is  M-step  reachable.  Then  from 
the  augmentation  results  of  the  single-input  case  in  Chapter  Three, 
there  exist  an  (M-n)  x n matrix  D,  an  (M-n)  x (M-n)  matrix  B and  an 
(M-n)  X 1 vector  e such  that  the  system  (F,g)  is  M-step  reachable  where 


A 

F = 

■ F-GLg 

0 ' 

> 

II 

< 

Gug 

D 

B 

e 

In  this  case,  the  closed-loop  feedback  system  for  (F,G)  has  the 
structure  in  Fig.  5.1,  where  [Aj  is  the  feedback  gain  matrix  for 

the  single-input  augmented  system  (F,g)  and  Ug(k)  = -Jli(k)x(k) 
-Jl2(k)w(k) . 
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E 

0) 


CO 

oo 

c 

Q) 

> 


> 


Fig.  5.1  Augmentation  for  M-cyclizable  systems 
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We  have  so  far  considered  cyclizabi lity  of  N-step  reachable  time- 
varying  systems.  But,  we  can  not  use  this  result  to  construct  a 
feedback  system  for  arbitrary  N-step  reachable  systems  since  it  is 
unknown  whether  or  not  N-step  reachability  always  implies 
M-cyclizability  with  M > N > n. 

We  now  consider  cycl izabi 1 ity  of  the  augmented  system  constructed 
in  Chapter  Three.  In  the  following,  we  shall  show  that  the  augmented 
system  constructed  in  Chapter  Three  is  always  cyclizable. 

(5.4)  THEOREM.  Suppose  that  an  m-input  n-dimensional  system  (F,G)  is 
N-step  reachable  with  N > n and  that  for  some  integer  q where 

A A 

N < q < mN,  (F,G)  is  the  q-dimensional  augmented  system  constructed  in 
the  proof  of  Theorem  3.9.  Then  (F,G)  is  q-cycl izable.  That  is,  there 
exist  a matrix  Lq  e and  a vector  ug  e A"’  such  that  (F-GLq,  Gu^)  is 
q-step  reachable. 


PRO^.  For  the  simplicity  of  the  proof,  let  us  choose  g^  as  a 
generator.  Then  uq  = [1  0 ...  0]'  and  Gug  = g^^ . From  the  generalized 
controllability  indices  determined  in  the  proof  of  Theorem  3.10,  let 

« - r.,  \ 

“max  ■ 1 < i < 


In  the  following,  we  let  ng  = 0. 
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We  now  define  integers  t and  as  follows: 

Initialization:  t=0,  i=l,a=0. 

— ..^P  ^ ^ J’t  “ ^ ^ 1 ) » increase 

£=1 

t by  1,  and  go  to  Step  2. 

If  n^j  < i,  then  go  to  Step  2. 

^ • If  a < m,  then  increase  a by  1 and  go  to  Step  1. 

If  a = m and  if  i < then  increase  i by  1 and  go 

to  Step  1. 

If  a = PI  and  if  i = a,^g^,  then  we  end  this  algorithm. 


From  the  integers  1 < t < q,  1 < < q,  we  let  P be  a q x q matrix 

'f  h 

whose  t'-^-column  vector  is  a q x 1 unit  vector  having  1 at 
position.  The  matrix  P is  called  a permutation  matrix. 

From  the  coefficients  r]''  e A of  (3.12),  we  define  the  elements 
e A as  fol lows: 


a 


z 

n. 


^ j+1  ’ 


1 

0 


, 1 < t < j 

z = j+1 

j+2  < t < m 


for  Z = 1,2,  . . 


fn  ; j = 1,2,  . . . , m-2 


m-1 


£ = 1,2, 


m-1 


1 


, Jl  = m. 


Now  define  the  vectors  u^-j  £ A'^  f or  i = 1 , 2 , . . . , m-1 , j=l,2, 
. . . , n-j  as  follows  ; 


< 


li  1 


n-j 


, Til 
Hi-i 


when  j = n^ 


otherwi se 


For  i = 1,  2,  . . . , tn-1,  j = 1,  2,  . . . , n,  let  B = ( J]  n.  , ) + 

1=1 

Define  Ug  = for  1 < i < m-1 , 1 < j < n^-  and  Ug  = 0 for  0^  + 02  + 
• * • ^m-1  + 1 < 3 < q-1. 

Define  Cq»  . . . , as 


Co  = Guq 

^i  = + 6u. , i = 1,  2,  . . . , q-1. 


Then  it  can  be  shown  that 
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where 


Cq  - [g^  . 


n ^ " 1 A 

s.  92 


and  Uq  e is  the  upper  triangular  matrix  whose  diagonal  elements 

I's. 


are 


Since  P,  Cq,  Uq  have  rank  q at  all  times, 


rank  [Cq  ...  Cq_i]  = q at  all  times. 


Let  Lg  - -[uj  U2  . . . Uq_2  0]a[Cg  Then  the 

reachability  matrix  of  the  system  (F-GLq,Guq)  is  [Cg  . . . ?q_i]. 

A A 

Therefore,  (F-GLg,  Gug)  is  q-step  reachable  and  thus  the  theorem  is 
proved. 


Again  let  (F,G)  be  an  N-step  reachable  system  with  reachability 
matrix  Cf^(k),  and  consider  the  array  of  vectors  shown  below  : 


91 

92 

• • • 9m 

Sf9i 

• 

%92 

• 

• • • Spgp^ 

• 

• 

5f  9i 

• 

• 

Sf  gi 

• 

• 

...  ^F  9[ 

By  searching  the  columns  of  this  array  from  top  to  bottom,  we  select 


those  vectors  that  are  not  equal  to  an  A-1 inear  combination  of 
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previously  selected  vectors.  Vectors  that  are  A-1 inear  combinations  of 
previously  selected  vectors  are  discarded.  In  this  selection,  we  assume 
that  there  exist  an  integer  t (where  1 < t < m)  and  integers  ni.n?, 

. . . , n^  (where  1 < n < N f or  i = 1 .2 t)  such  that 


(5.5) 


I 

j=l 


"j 

I 

1=0 


(S^gj) 


1,  2, 


t 


are  elements  in  A.  and  that  the  vectors  Oi . Spgi  , . . . , 

3.1  ? 92 > • ♦ • > » • » « » Sp  q^■  generate  A'^.  Here,  g^^,  g2» 

. . . , g^  are  all  included  in  our  selection.  If  this  is  not  the  case, 
then  we  can  rearrange  the  above  array  by  reordering  those  columns. 
Therefore,  we  can  wok  with  the  reordered  array.  Letting  q = n^  + n2  + 

. . . + n^,  we  have  N < q < mN. 

If  the  system  (F,G)  is  time-invariant,  there  always  exist  integers 
t,  nj,  n2,  . . . , n^  such  that  (5.5)  holds  (see  KAILATH  [1980]). 
However,  in  the  time-varying  case,  there  may  not  exist  Sp'^g^-  for  all  1 < 
i < m and  for  all  1 < j < N such  that 


S 


j 

F 


9i 


i=l  j-1 


I I (Sp^g.) 

r=l  i=0  ^ ^ 


e A 


Therefore,  we  are  considering  a very  restricted  case  of  Theorem  5.4. 
However,  this  restriction  gives  a very  simple  constructive  proof  of 
cycl izabi 1 ity . 

We  now  consider  the  augmentation  for  the  system  (F,G)  which  is 


assumed  as  above. 
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Choose  q vectors  1 < j < t,  0 < i < rij-l,  where  e^-j  are 

(q-n)  X 1 vectors,  such  that 


rank 


r.  c 

9l  • • • 9l  92  • • . 9t  • • • Sp  ^ ■ 

®01  • • • ®{ni-l)l  ®02  • • • '^Ot  • • • ®(n^-l)t 

at  all  times. 


q 


This  is  always  possible  from  Proposition  3.4.  Let  Cq(k)  denote  the 
q X q matrix 


(5.6)  Cq(k)  = 


9l  • • • 9l  92  . . . 9t  . . . % ^ 9t 

_®01  • • • e(ni-l)l  eo2  . . . eot  . . . 


From  (5.5)  and  (5.6),  calculate  e^.i,  i = 1,2,  . . . ,t  , such  that 
(5.7) 


-nil 


i , 

I I e . . 

j=l  1=0 


Define  E(k)  = [eoi(k)  eQ2(k)  . . . eot(k)  * . . . *]  and 


[D(k)  B(k)]  = X(k)(aC-M(k) 


where  *'s  are  arbitrary  vectors  and  X(k)  = [en  . . . &i2 

• • • ®n^t^* 
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Let  F = 

U_ 

1 

1 

o 

and  G = 

“ — 

G 

1 

o 

CO 

E 

Then  by  the  same  argument  in 


the  proof  of  theorem  3.10,  we  see  that  the  q-dimensional  augmented 

A A 

system  (F,G)  is  q-step  reachable  at  all  times.  In  the  following 
theorem,  we  shall  show  that  the  augmented  system  constructed  above  is 
q-cycl izable. 


(5.8)  THEOREM.  Let  (F,G)  be  an  N-step  reachable  system  having  the 
property  (5.5).  Then  for  some  integer  q where  N < q < mN  the  augmented 

A A 

system  (F,  G)  constructed  above  is  q~cyclizable,  i.e.  there  exist  a 
matrix  Lg  s and  a vector  ug  c such  that  (F  - GLQ,Gug)  is  q-step 
reachable. 


PR00F_.  Without  loss  of  generality,  we  choose  g^  as  a generator.  Then 
Uq  = [1  0 . . . 0] ' . Oef i ne 

A 

Co  = Gug 

= F(ct5^-_j)  + Gu^-,  i = 1,  2,  . . . , q-1 

where 

^i  “ [0  ...  1 ...  0 ] ‘ , when  i = n^ , n2  , . . . , n^ 

(i-n  +1)^^ 

' G , otherwise. 
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Then  we  have 

A 

Co  = Guq  = 


n^-l. 


9l 


F(aC 


n^.l)  + 92  = 92  + S 


^q-n^+1  9^  t ‘ * 

Ht-1  . 

^q-1  9^  + . . . 


If  we  rewrite  the  above  vectors  in  a matrix  format,  then  we  obtain 


[?0  Cl  . . . Cq_i]  = Cq  Uq 

where  Cq  is  given  by  (5.6)  and  Uq  is  an  upper  triangular  matrix  whose 
diagonal  elements  are  I's.  Since  Cq,  Uq  are  invertible,  [^q  Ci  . . . 
?q_l]  is  also  invertible.  Define 


Lq  - - [ui  U2  ...  Uq_i  0 laC^o  ?1  • • • ^q-l]'^  • 

Then  the  reachability  matrix  of  the  system  (F-GLq,Guq)  is  [Cq  Ci  . . . 
Cq_i]  and  thus  the  theorem  is  proved.  [] 
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We  now  construct  a feedback  gain  matrix  for  the  augmented  system 

A A 

(F,G)  by  the  method  discussed  in  this  chapter.  Since  the  single-input 

A A A 

system  (F-GLq,Guq)  constructed  in  Theorem  5.4  and  Theorem  5.8  is  q-step 
reachable,  from  (2.8)  there  exists  a coordinate  transformation  T(k)  such 
that 


F(k)  = T(k+l)(F(k)  - G(k)L  (k))T-l(k) 

g(k)  = T(k+l)Guo 

are  in  the  control  canonical  form.  Let  us  define  ag,  aj,  . . . , e 
A as 

(5.9)  [ag  ai  . . aq_i]'  = -[Cg  . . . ?q_i]-l  (S  . . (Gu.))  . 

(F-GLq)  ° 

We  pick  any  monic  polynomial  x(z)  of  degree  q as  a characteristic 
polynomial  of  the  closed-loop  system  for  (F-GLq,Guq),  where  x(z)  has  the 
form 

x(z)  = z9  + I z^'b.  , b.  e R . 

i=0  ^ ^ 

Then  t(k)  - -[bg  - ag  . . . bq_j  - aq_i]T(k)  is  the  feedback  gain  matrix 
for  the  single-input  system  (F-GLq,Guq)  which  makes  the  system  have  a 
characteristic  polynomial  x(z).  Therefore,  if  we  let 
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(5.11)  L(k)  = Lo(k)  - uoA(k)  , 

then  L(k)  is  the  feedback  gain  matrix  for  the  augmented  system  (F,G), 
and  the  closed-loop  system  will  have  the  characteristic  polynomial  x(z). 
We  illustrate  the  above  results  by  the  following  example. 


(5,12)  EXAMPLE . Consider  the  system  (F,G)  such  that 


F(k)  = 

1 0 

and 

G(k)  = 

’ 1 o' 

o 

Cl. 

1 

ev(k)  od(k) 

Then  the  reachability  matrix  is 


C2(k) 


1 0 1 
ev(k)  od(k)  0 


0 

ev(k) 


and  this  system  is  2-step  reachable.  We  select  g^ , g2»  Spg^  from  C2(k) 
as  the  vectors  that  are  not  equal  to  an  A-1 inear  combination  of 
previously  selected  vectors.  Then  n^  = 2,  n2  = 1 and  q = 3.  Now 
determine  eg^,  eQ2,  such  that 


C3(k)  = 


91 

eOl 


has  rank  3 at  all  times. 


92 

Sf9i 

®02 

^11 

= ev(k),  e^i  = 0,  we  have  rank  C3(k)  = 3 at 


Letting  egi  = -od(k),  eQ2 
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all  times.  From  Sp^g^  = -93  + Spg^  and  Sp92  = gj  - Spg^,  calculate  e2i 
and  e-^2>  Then  we  have 

621  ^ -od(k) 

612  = -6v(k) 


Define  E(k)  = [egi  eQ2]  = [-od(k)  ev(k)]  and 


[D(k)  B(k)]  = [eii  e2i  6i2](aC3)"^(k) 

= [0  -od(k)  -ev(k)] 


= [-ev(k)  0 -1  ]. 


0 od ( k ) 

0 ev(k) 

1 -od(k) 


-ev(k) 

od(k) 

ev(k) 


Then  the  augmented  system  (F,G)  has  the  form 


1 

0 

0 ■ 

and 

A 

G = 

1 

0 ■ 

-od(k) 

1 

0 

ev(k) 

od(k) 

. -6v(k) 

0 

-1  _ 

-od(k) 

ev(k) _ 

Now,  let  us  construct  a feedback  control  system  for  the  augmented  system 
(F,G).  Let 


UQ 


1 

0 


75 


Then  Cq  = Guq  = = 


1 

ev(k) 
-od(k) 


Take  uj 


-1 

and  U2  = 

1 

0 

1 

_ 0_ 

_2_ 

A 

A 

Cl  = F(crg^)  -g^ 


= F(aCj)  + 92  ^ 


Then  we  have 


0 

-ev(k) 
od  ( k ) 
0 

od(k) 

ev(k) 


Since  [^g  ?i  ^2]  is  invertible,  we  can  define 


Lo  = - [ui  U2  0 ]aC5g  ^2^"^ 


-100' 

’ 1 0 o' 

0 2 0 

1 -od(k)  ev(k) 

- 0 ev(k)  od(k)_ 

= 1 0 0 

.-2  2od(k)  -2ev(k)_ 

From  (5.9),  we  have  ag  = 0 , ai  = 1,  a2  = 0.  Therefore,  the  coordinate 
transformation  matrix  T(k)  will  be 
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T(k)  = 

1 

-od(k) 

ev(k ) 

1 

-ev(k) 

ev(k ) 

0 

ev(k) 

od(k) 

Pick  x(z)  = z^  - zO.25.  Then  from 

(5.11)  we 

have 

L(k)  = Lq  + ugChg-ag  bj-aj  b2-a2]T(k) 

= 0.75  0.25od(k)  -0.25ev(k) 

. -1  od(k)  -ev(k) 


If  we  construct  the  feedback  gain  matrix  from  the  generalized  control 
canonical  form  described  in  Chapter  Three  and  Chapter  Four,  then  we  have 


L(k) 


1 + 1.75(od(k)  - ev(k))  -1.75od(k)  -1.75ev(k) 

_od(k)  - ev(k)  -od(k)  -ev(k) 


The  simulation 
and  the  initial 


results  are  given  in  Fi 
conditions  are  x(0)  = 


g.  5.2 


2 

1 


and  Fig.  5.3 
, w(0)  = 0.2 


when  v(k) 


=0 


It  is  interesting  to  compare  the  responses  of  the  closed-loop  systems 
constructed  by  those  two  methods. 
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a)  States  x(k)  = [x^(k)  X2(k)]'  and  w(k) 


1 

2 


b)  Input  u(k)  = [u^(k)  U2(k)]' 


Fig.  5.2  Simulation  results  of  Example  5.12  (Canonical  form 
method:  v(k)  = 0,  x^(0)  = 2,  X2(0)  = 1,  w(0)  = 0.2) 
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a)  States  x(k)  and  w(k) 


1 

2 


b)  Input  u(k) 


Figure  5.2  Simulation  results  of  Example  5.12  (Direct  method: 
v(k)  = 0,  Xj(0)  = 2,  X2(0)  = 1,  w(0)  = 0.2) 


CHAPTER  SIX 

TIME-VARYING  RESULTANT  MATRICES 

Th0  second  part  of  this  dissertation  is  concerned  with  developing 
the  theory  of  the  resultant  for  linear  time-varying  systems.  The 
resultant  is  a classical  tool  used  in  examining  the  relative  primeness 
of  a prescribed  pair  of  polynomials.  In  the  time-invariant  case,  the 
resultant  has  a wide  range  of  applications  in  system  theory  especially 
when  the  system  is  given  by  a polynomial  matri x-f racti on  description. 
For  example,  by  obtaining  greatest  common  divisors  of  polynomials  and 
polynomial  matrices  via  the  resultant,  we  can  determine  uncontrollable 
and  unobservable  modes  of  given  systems. 

However,  little  attention  has  been  paid  to  the  development  of  the 
resultant  in  the  time-varying  case  since  a polynomial  matrix-fraction 
description  for  time-varying  systems  was  not  available  until  KAMEN 
et  al . [1985]  and  POOLLA  [1984]  successfully  developed  a time-varying 
version  of  the  polynomial  matrix-f raction  approach  in  the  framework  of 
skew  rings  of  polynomials,  formal  power  series,  and  formal  Laurent 
series.  Moreover,  their  results  closely  resemble  the  time-invariant 
case.  Using  their  results,  we  define  the  resultant  for  time-varying 
systems  and  exploit  several  applications  to  system  theory  in  the  second 
part  of  this  dissertation. 
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In  this  chapter,  we  define  a time-varying  version  of  Sylvester's 
resultant  matrix  for  polynomials  and  polynomial  matrices.  It  will  be 
shown  that  reachability  of  the  observable  realization  (4.12a,b)  and 
(4.18a,b)  can  be  checked  by  the  resultant  matrix.  We  begin  with  the 
definition  of  the  resultant  matrix  for  two  polynomials  with  time-varying 
coefficients. 

Let  a pair  of  polynomials  with  time-varying  coefficients,  a(z)  and 
b(z),  be  given  by 

a(z)  = z"^  + z^“^  + . . . + zai(k)  + ag(k) 

(6.1) 

b(z)  = + . . . + zbi(k)  + bo(k)  . 

Here,  we  should  note  that  a(z)  and  b(z)  are  elements  of  right  skew 
polynomial  ring  A[z]  defined  in  Chapter  Two.  The  indeterminate  z has 
the  noncommutati ve  property  az  = z(aa),  a e A,  where  o is  the  right- 
shift  operator  such  that  (aa)(k)  = a(k-l).  Then  we  have  the  following 
time-varying  version  of  the  resultant  matrix. 

(^•^)  definition . The  time-varying  resultant  matrix  rp(a,b)  of  two 
polynomials  a(z)  and  b(z)  given  by  (6.1)  is  defined  as 


81 


rn(a,b) 


ao  0 

aj  aao 

• • 

^n-1 

1 oa^_i 

1 • 


a ag 


_n-2, 

^n-l 

1 


0 

abg 


'^n-l 
0 ab 

0 
0 


n-1 


^n-lk 
a Dr 


a'l-lb 


n-1 


where  rp(a,b)  e ^ (2n-l) 


[] 


The  above  definition  closely  resembles  Sylvester's  resultant  matrix 
in  the  time-invariant  case  (see  BARNETT  [1983]).  In  the  following,  we 
relate  the  definition  of  the  time-varying  resultant  matrix  to 
reachability  of  a given  system. 

Let  us  consider  the  input/output  difference  equation  given  by 

(6.3)  y(k+n)  + I ai(k)y(k+i)  = I b,- (k)u(k+i')  . 

i=0  i=o 

The  observable  realization  (F,g,h)  of  equation  (6.3)  is  given  by 

(4.18a,b).  Define  polynomials  related  to  the  input/output  difference 
equation  (6.3)  as  follows: 


82 


a(z)  = z"  + zn-lan_i(k-n+l)  + . . . + zai(k-l)  + ao(k) 

(6.4) 

b(z)  = z^"^b^_i(k-n+l)  + . . . + zbi(k-l)  + bo(k)  . 

Then  reachability  of  the  observable  realization  (F,g,h)  of  (6.3)  can  be 
examined  by  the  resultant  matrix  of  the  two  polynomials  a(z),  b(z)  given 
by  (6.4). 

(6.5)  THEOREM.  The  n-dimensional  observable  realization  (F,g,h)  given 
by  (4.18a,b)  is  n-step  reachable  at  all  times  if  and  only  if 

det  rp(a,b)  f 0 at  all  times 

where  a(z),  b(z)  are  given  by  (6.4)  and  det  denotes  the  determinant  of  a 
matrix. 


PROOF.  First,  partition  the  resultant  matrix  r^(a,b)  as  follows: 

n-1  n 


Tn(a,b) 


Postmultiply  A = 


^1 


n 

n-1 


on  both  sides  of  the  above  equation. 


where  is  an  (n-1)  x n matrix  satisfying  the  following: 


+ T4  = 0 
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Then  we  have 
(6.6)  r^(a,b)A  = 

Since  F3  is  an  invertible  (n-1)  x (n-1)  matrix,  there  always  exists  a 
matrix  Aj  such  that  r^A^  + = 0.  It  remains  to  show  that 


where  = [g  Spg  . . . Sp"^g].  This  is  done  by  comparing  columns  on 
both  sides  of  the  above  equation  with  A^  = -r3~^r4.  From  (6.6),  it  can 
be  shown  that 


FiAi  + 


det  r^(a,b)  = det  A det 

= (-l)"("-')/2det  C„  . 

Hence,  the  system  (F,g,h)  is  n-step  reachable  at  all  times  if  and  only 
if  det  r^(a,b)  0 at  all  times.  [] 

In  the  time-varying  case,  an  n-dimensional  system  can  be  N-step 
reachable  with  N > n.  Therefore,  we  need  to  generalize  Theorem  6.5  to 
the  N-step  reachable  case.  To  do  this,  we  modify  Definition  6.2  as 


fol 1 ows: 
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(6.7)  DEFINITION.  For  an  integer  N > n,  the  time-varying  resultant 
matrix  r^(a,b)  of  two  polynomials  a(z)  and  b(z)  given  by  (6.1)  is 
defined  as 


fN(3,b)  = 


ao  0 

ai  aag 


a aQ 


^n-l 

1 aa 
1 


n-1 


• • • 


1 


t>0 


0 

abo 


’n-1 
0 ab 


n-1 


a bg 


aN-lb 


n-1 


where  r[^(a,b)  e ^ (2N-1)^ 


[] 


When  M=n,  the  above  definition  is  exactly  same  as  Definition  6.2. 
But,  if  N > n , then  the  resultant  matrix  is  no  longer  a square  matrix. 

(6.8)  THEOREM.  The  n-dimensi onal  observable  realization  (F,g,h)  given 
by  (4.18  a,b)  is  N-step  reachable  with  N > n if  and  only  if 


rank  r|^(a,b)  = N + n - 1 at  all  times 
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where  a(z)  and  b(z)  are  given  by  (6.4). 


PROOF.  This  proof  is  similar  to  that  of  Theorem  6.5.  First,  we 
partition  the  resultant  matrix  r[^j(a,b)  as  follows: 


Postmultiply  A = 


rf^(a,b)  = 


n ^N-1 


N-1 
■ ^1 


N 


P4 


n 

N-1 


on  both  sides  of  the  above  equation, 


where  Aj^  - -r3“^r4.  Then  we  have 


TM(a,b)A 


+ 


P2 


0 


It  can  be  shown  that  FjA;^  ^2  ” 0|\j  by  comparing  columns  on  both  sides, 
where 


On  =1^9  Spg  . . . Sp^“^g] 

Since  rank  F3  = N-1  and  rank  A = 2N-1  at  all  times,  rank  r^A^  + T2  = n 
at  all  times  if  and  only  if  rank  = N + n - 1 at  al  1 times. 

Therefore,  the  theorem  is  proved.  [] 
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Now  let  us  consider  the  MIMO  case.  First,  we  define  a time-varying 
version  of  ROWE's  [1972]  resultant  matrix.  Later  we  will  find  a 
generalized  resultant  matrix  which  is  analogous  to  that  of  ANDERSON  and 
JURY  [1976]. 

Let  polynomial  matrices  with  time-varying  coefficients,  Q(z)  and 
R(z),  be  given  by 


where  are  p x p and  R^  are  p x m time-varying  matrices. 

(6.10)  DEFINITION.  The  time-varying  resultant  matrix  rp(Q,R)  of  two 
polynomial  matrices  Q(z)  and  R(z)  given  by  (6.9)  is  defined  as 


Q(z)  = z^I  + z^"^Qn_i(k)  + . . . + zQi(k)  +Qo(k) 


(6.9) 


R(z)  = z^"k^_;^(k)  + . . . + zRi(k)  + Rg(k) 


Qo  0 

“ Ql  ®Qo 


RO  0 

Rl  oRq 


0 


0 


0 
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where  r„(Q,R)  t A<"P  - DP  » Knp-Dp  + r.pm>.  [] 

Suppose  that  the  MIMO  linear  time-varying  discrete-time  system  is 
given  by  the  input/output  difference  equation 

Pl“l  P2“l 

(6.11)  y(k+pi)  + I Qi(k)y(k+i)  = I Ri(k)u(k+i) 

i=0  i=0 

where  pj  > P2»  Qi  e AP^P,  and  e AP^’’’.  Then  the  n-dimensi onal 
observable  realization  (F,G,H)  of  the  above  equation  is  given  by 
(4.12a,b),  where  n = pjp  . Define  polynomial  matrices  related  to  (6.11) 
as  follows  : 


Q(z)  = z*^^I  + z'^^  Qpj-l('^"Pl‘'’l)  + . . . + zQi(k-l)  + Qo(k) 

(6.12) 

R(z)  = z^^  f^p2-l (*^~P2''‘^ ) + . . . + zRi(k-l)  + Rg(k) 

Then  we  have  the  following  theorem. 

(6.13)  THEOREM.  The  n-dimensional  observable  realization  (F,G,H)  given 
by  (4.12a,b)  with  n = pjp  is  n-step  reachable  at  all  times  if  and  only 
if 

rank  rp^(Q,R)  = (p^p  + Pi  - l)p  at  all  times  , 


□ 


where  Q(z)  and  R(z)  are  given  by  (6.12). 
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We  omit  the  proof  of  the  above  theorem  since  it  is  almost  identical 
to  that  of  Theorem  6.5.  As  in  the  SISO  case,  we  can  define  the  N-step 
time-varying  resultant  matrix  in  the  MIMO  case. 


(6.14)  DEFINITION.  For  an  integer  N > p the  time-varying  resultant 
matrix  F|^(Q,R)  of  two  polynomial  matrices  Q(z)  and  R(z)  given  by  (6.9) 
is  defined  as 


Fn(Q.R)  = 


Qo  0 

Ql  cjQo 


RO  0 
Rl  ctRq 


Qn-1 

I crQ^.i 
I 

0 


Rn-l 

0 aR^.i 


I 0 


where  r^(Q,R)  e x {(Np-l)p  + Npm}. 


[] 


(6.15)  THEOREM.  The  n-dimensi onal  observable  realization  (F,G,H)  given 
by  (4.12a,b)  with  N > n = pjp  is  N-step  reachable  at  all  times  if  and 
only  if 
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rank  r[>^(Q,R)  = (Np  + n - l)p  at  all  times  , 


where  Q(z)  and  R(z)  are  given  by  (6.12). 


[] 


The  shortcoming  of  the  time-varying  ROWE's  resultant  is  that  the 
size  of  the  matrix  is  too  big  to  handle.  Therefore,  we  need  to  reduce 
the  size  of  the  resultant  matrix.  In  the  time-invariant  case,  ANDERSON 
and  JURY  [1976]  have  defined  a generalized  Sylvester's  resultant 
matrix.  We  will  obtain  the  similar  generalized  resultant  matrix  in  the 
time-varying  case. 

(6.16)  DEFINITION . For  each  integer  i > 0 , a generalized  time-varying 
resultant  matrix  of  two  polynomial  matrices  Q(z)  and  R(z)  given  by  (6.9) 
is  defined  as 


0 


0 


ri(Q,R)  = 


0 


0 


0 


0 


% ^ 


0 
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where  rMQ,R)  0 + n - l)p  x {(i  - l)p  + im)^ 

Using  the  above  definition,  we  get  a necessary  and  sufficient 
condition  for  n-step  reachability  of  the  given  system. 

(6,17)  THEOREM.  The  n-dimensional  observable  realization  (F,G,H)  given 
by  (4.12a,b)  with  n = p^p  is  n-step  reachable  at  all  times  if  and  only 
if  there  exists  an  integer  a < np  such  that 

(i)  r“(Q,R)  has  full  rank  at  all  times 

(ii)  rank  r“+l(Q,R)  - rank  r«(Q,R)  = p at  all  times. 

The  proof  of  the  above  theorem  can  be  done  by  a slight  modification 
of  that  of  ANDERSON  and  JURY  [1976].  Reachability  of  the  system  (F,G,H) 
can  be  checked  by  examining  the  ranks  of  generalized  resultant  matrices 
of  growing  size.  Therefore,  the  computation  may  be  reduced  when  the 
generalized  resultant  matrix  is  used  for  checking  reachability  of  the 
system,  we  should  note  that  if  a > n,  then  Definition  6.16  is  exactly 
same  as  ROWE's  resultant. 


CHAPTER  SEVEN 

APPLICATIONS  OF  TIME-VARYING  RESULTANT  MATRICES 

Having  defined  the  discrete-time  time-varying  version  of  the 
resultant  matrix  in  the  previous  chapter,  we  now  continue  the 
development  of  the  resultant  theory  in  the  time-varying  case.  In  the 
following,  we  shall  consider  coprimeness  of  polynomials  and  polynomial 
matrices  with  time-varying  coefficients,  a time-varying  version  of 
Bezout  polynomial  factorization,  and  compensator  design  of  time-varying 
systems  using  the  resultant  matrix. 

Let  a pair  of  polynomials  with  time-varying  coefficients,  a(z)  and 
b(z),  be  given  by 

a(z)  = + z^"^  ^n-l^^)  + . . . + zai(k)  + agCk) 

(7.1) 

b(z)  = z"^  ^ti^_^(k)  + . . . + zbi(k)  + bg(k) 

If  there  is  a polynomial  c(z)  e A[z]  such  that  deg  c(z)  < n and  for  some 

polynomials  a(z)  e A[z]  and  b(z)  e A[z]  , a(z)  = c(z)a(z)  and  b(z)  = 
c(z)b(z),  then  the  polynomial  c(z)  is  said  to  be  a left  common  factor  of 

two  polynomials  a(z)  and  b(z).  If  a(z)  and  b(z)  do  not  have  a left 

common  factor,  a(z)  and  b(z)  are  said  to  be  left  coprime  (or  relatively 
left  prime).  Then  we  have  the  following  result. 
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(7.2)  theorem.  Two  polynomials  with  time-varying  coefficients,  a(z)  and 
b(z),  given  by  (7.1)  have  a left  common  factor  if  and  only  if 

det  r^(a,b)  = 0 at  all  times 

where  rp(a,b)  is  the  resultant  matrix  of  a(z)  and  b(z)  defined  by  (6.2). 

PRO^.  Suppose  that  the  polynomials  a(z)  and  b(z)  given  by  (7.1)  have  a 
left  common  factor  c(z)  of  degree  0 < r < n-1.  Then  c(z)  must  be  monic 
since  a(z)  is  monic.  Let 

a(z)  = c(z)a(z)  , b(z)  = c(z)b(z) 
where  a(z),  b(z),  c(z)  are  given  by 

a(z)  = z"---  + z"-'--laVr.i(k)  t . . . + ao(k) 
b(z)  = z"-'-lb„.,.i(k)  + . . . + bo(k) 
c(z)  = z''  + z'-lc^.iCk)  + . . . t co(k)  . 


Then  we  have 
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rn(a,b)  = 


Cq  0 
Cl  acQ 


3q  ° 

ai  aag 


Cp_l 

1 acp.i 


^n-r-1 

1 • . 


n-r- 


1 


By  Sylvester's  inequality,  rank  r^(a,b)  < (2n  - r -1)  at  all 
Therefore,  det  r^(a,b)  = 0 at  all  times. 

Conversely,  suppose  that  det  rn(a,b)  = 0 at  all  times, 
always  exists  (2n-l)-element  column  vector  x 0 such  that 


rp(a,b)x  = 0 at  all  times. 

Let  X - [cq  Cl  . . . Cp_2  dQ  di  . . . dp_i]'  and  define 

c(z)  = Cn_2  + . . . + zci  + Cq 
d(z)  = z^-1  d^.i  + . . . + zdi  + do  . 

Then  we  have 


0 

abQ 

1 


times . 

Then  there 


a(z)c(z)  + b(z)d(z)  = 0 
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Since  a(z)  is  monic,  it  is  invertible  over  A((z“M).  Hence  we  get 


c(z)  = -a“^ (z)b(z)d(z) 

The  polynomial  b(z)d(z)  must  have  a(z)  as  a left  factor  since  c(z)  is  a 
polynomial.  However,  degrees  of  b(z)  and  d(z)  are  n-1  and  n-1, 
respectively.  Therefore,  a(z)  and  b(z)  must  have  a left  common  factor 
of  degree  at  least  1.  [] 

Combining  the  above  result  and  Theorem  6.5,  we  see  that  the 
observable  realization  (F,g,h)  of  given  polynomials  a(z)  and  b(z)  is  not 
reachable  if  a(z)  and  b(z)  have  a left  common  factor.  But,  even  if 
(F,g,h)  is  not  reachable,  a(z)  and  b(z)  may  not  have  a common  factor. 

Ilchmann  et  al . [1984]  have  defined  a notion  of  coprimeness  of  two 
time-varying  polynomials  and  polynomial  matrices  in  terms  of  the  time- 
varying  version  of  Bezout  polynomial  factorization.  However,  this 
definition  of  coprimeness  is  not  related  to  a notion  of  the  resultant. 

In  the  following,  we  derive  the  time-varying  version  of  Bezout  identity 
from  the  reslutant  matrix,  and  we  relate  this  result  to  coprimeness  of 
time-varying  polynomials. 

(7.3)  THEOREM.  Let  the  polynomials  a(z)  and  b(z)  be  given  by  (7.1). 

Then  there  exist  polynomials  c(z)  and  d(z)  such  that 


(7.4) 


a(z)c(z)  + b(z)d(z)  = 1 
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where  deg  c(z)  < n-2,  deg  d(z)  < n-1  if  and  only  if 


det  r^(a,b)  ^0  at  all  times. 

£R00F_.  Let  c(z)  and  d(z)  be  denoted  by 

c(z)  = + . . . + ZCj  + Cq 

d(z)  = z'^"^d^_;^  + . . . + zdj  + dg 

If  we  rewrite  (7.4)  in  a matrix  format,  then  we  obtain 

Ccq  . . . Cf,_2  dg  . . . dn_i]  r^'Ca.b)  = [1  0 ...  0]  . 

The  above  equation  has  a solution  if  and  only  if  det  r^(a,b)  0 at  all 

times,  and  thus  the  theorem  follows.  [] 

(7.5)  COROLLARY . Let  the  polynomials  a(z),  b(z)  be  given  by  (7.1). 

Then  for  some  integer  N > n there  exist  the  polynomials  c(z)  and  d(z) 
such  that 

a(z)c(z)  + b(z)d(z)  = 1 


where  deg  c(z)  < N-2,  deg  d(z)  < N-1  if  and  only  if  rank  r|y|(a,b)  = N + n 
- 1 at  all  times.  n 
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From  the  above  results,  we  see  that  the  existence  of  a Bezout 
identity  is  necessary  but  not  sufficient  for  coprimeness  of  two 
polynomials. 

(7.7)  SAMPLE . Consider  the  polynomials  a(z)  and  b(z)  given  by 

a(z)  = z^  + zk  + k 
b(z)  = z 


Then  we  have  the  resultant  matrix  r2(a,b)  given  by 


F2(a,b) 


k 0 

k 1 

1 0 


0 

0 

1 


Since  det  r2(a,b)  t 0 for  all  integers  k except  k = 0,  a(z)  and  b(z)  are 
left  coprime.  However,  for  any  integer  N > 2,  rank  r,^(a,b)  N + 2 - 1 

at  all  times.  Therefore,  a(z)  and  b(z)  does  not  satisfy  a Bezout 
identity.  |-^ 

Having  studied  the  results  of  coprimeness  and  Bezout  identity  in 
the  SISO  case,  we  now  consider  the  MIMO  case. 

Let  polynomial  matrices  with  time- varying  coefficients,  Q(z)  and 
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R(z),  be  given  by 

Q(z)  = z^l  + + . . . + zQi  + Oq 

(7.8) 

R(z)  = z^'^R^.i  + . . . + zRi  + Rg 

where  e AP^P  and  R^-  e AP^*^.  Here,  we  assume  that  p < m. 

(7.9)  THEOREM.  The  polynomial  matrices  Q(z)  and  R(z)  given  by  (7.8) 
have  a left  common  factor  if  and  only  if  for  any  integer  i > 0 

rank  rMO.R)  < (i+n-l)p  at  all  times  [] 

(7.10)  THEOREM.  Let  the  polynomial  matrices  Q(z)  and  R(z)  be  given  by 
(7.8).  Then  for  some  integer  i > 2 there  exist  polynomial  matrices  C(z) 
and  D(z)  such  that 

Q(z)C(z)  + R(z)D(z)  = I 

where  deg  C(z)  < i-2,  deg  D(z)  < i-1  if  and  only  if 

rank  r^Q.R)  = (i  + n - l)p  at  all  times  [] 


As  we  have  mentioned  in  the  SISO  case,  the  existence  of  Bezout 
identity  implies  the  coprimeness,  but  the  converse  is  not  always  true. 
This  situation  may  occur  since  even  if  the  polynomial  matrices  Q(z)  and 
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R(z)  are  left  coprime,  the  generalized  resultant  matrix  fMQ.R)  for  any 
integer  i may  not  have  full  rank  at  all  times. 

From  the  above  results  and  Theorem  6.17,  we  see  that  the 
n-dimensional  observable  realization  (F,G,H)  of  two  polynomial  matrices 
Q(z),R(z)  given  by  (6.12)  is  N-step  reachable  with  N > n if  and  only  if 
there  exist  the  polynomial  matrices  C(z)  and  D(z)  such  that 

Q(z)C(z)  + R(z)D(z)  = I 

where  deg  C(z)  < N-2,  deg  D(z)  < N-1.  This  result  corresponds  to  that 
of  POOLLA  [1984]. 

We  now  apply  the  resultant  matrix  to  the  design  of  output  feedback 
controllers  and  obtain  a result  on  a type  of  "assignability". 

Again  suppose  that  the  system  is  given  by  the  input/output 
difference  equation 

Pi"!  P2“l 

(7.11)  y(k+pi)  + I Q^(k)y(k+i)  = I R-(k)u(k+i) 

i=0  i=o 

where  p > P2,  e AP^P,  R^.  e aP^'^.  The  polynomial  matrices  related  to 
this  system  is  defined  by 

Q(z)  = z^^I  + z^l  ^p^_i(k-pi+l)  + . . . + zQi(k-l)  + Qg(k) 

pp-1 

^ Rp2-l(l<*P2'^l)  + . . . + zRi(k-l)  + RQ(k) 


(7.12) 
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Since  Q(z)  Is  invertible  over  A((z’M),  we  can  define  W(z)  as 

(7.13)  w(z)  = Q-l(z)R(z)  . 

Then  W(z)  is  said  to  be  the  transfer-function  matrix  of  the  input/output 
difference  equation  (7.11),  and  (7.13)  is  called  the  left  polynomial 
matrix-fraction  representation  of  (7.11)  (see  KAMEN  et  al . [1985],  and 
POOLLA  [1984]).  If  for  some  matrices  Q(z)  and  R(z),  a transfer-function 
matrix  W(z)  is  given  by 

W(z)  = R(z)Q-l(z)  , 

then  it  is  called  the  right  polynomial  matrix-fraction  representation. 

If  we  let  U(z),Y(z)  be  the  polynomial  representation  of  input  and  output 
sequences,  then  we  have 

(7.14)  Y(z)  = W(z)U(z) 

Consider  the  feedback  control  system  shown  in  Fig  7.1.  Here,  W^(z) 
is  the  transfer  function  of  the  compensator  to  be  designed,  which  has 
the  form 


W(.(z)  = C(z)D"^(z) 
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Fig  7.1  Feedback  control  system  structure. 


In  the  following,  we  design  a compensator  Wj-(z)  such  that  the 
closed-loop  system  has  preassigned  dynamics. 

Define 

(7.15)  E(z)  = W(-(2)Y(z) 

and  let  C(z)  be  a vector  such  that 

(7 .16  a)  D(z)5(z)  = Y(z) 

(7.16  b)  E(z)  = C(z)5(z) 

The  polynomial  vector  ?(z)  is  called  the  partial  state  of  the  system. 
From  (7.13)  and  (7.14),  we  have 

(7.17)  Q(z)Y(z)  = R(z)U(z). 

Let  V(z)  be  a reference  input.  Then  we  obtain 


V(z)  = U(z)  + E(z) 
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Premultiplying  R(z)  on  both  sides  of  the  above  equation  and  from  (7.16), 
(7.17),  we  obtain 

R(z)V(z)  = R(z)U(z)  + R(z)E(z) 

= Q(z)Y(z)  + R(z)C(z)C(z) 

= (Q(z)D(z)  + R(z)C(z))Uz). 

Inserting  the  above  equation  into  (7.16a),  we  can  get 

(7-18)  Y(z)  = D(z)(Q(z)D(z)  + R(z)C(z))-1r(z)V(z)  . 

Let  us  denote  n(z)  = Q(z)D(Z)  + R(z)C(z).  Then  we  see  that  the  closed- 
loop  system  dynamics  are  determined  by  the  polynomial  matrix  n(z). 

We  now  consider  assigning  the  coefficients  of  the  polynomial  matrix 
n(z)  by  choosing  C(z)  and  D(z). 

Suppose  that  the  n-dimensional  observable  realization  (F,G,H)  of 
the  polynomials  given  by  (7.12)  with  n = p^p  is  N-step  reachable  at  all 
times  with  N > n.  Then  from  Theorem  6.17,  there  exists  an  integer  a > 0 
such  that 

(i)  r“(Q,R)  has  full  rank  at  all  times 


(ii)  rank  r«‘^^(Q,R)  - rank  r«(Q,R)  = p at  all  times. 
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Let  the  polynomial  matrix  n(z)  be  given  by 

n(z)  = z'^l  + + . . . + zn^  + Hq 

where  v = (a  + n - 1 ) . 

Since  r®(Q,R)  has  full  rank  at  all  times,  there  always  exist 
matrices  Dg,  Dj,  . . . , D„_2  e Cg,  Ci,  . . . , s such 

that 


(7.19) 


Let  us  define 


r^(Q,R) 

■ d«-2  ■ 

• 

- 

’ ^v-1  ■ 

• 

• 

Do 

• 

• 

• 

Co 

• 

• 

.^a-1. 

. "0 

C(z)  = + . . . + zCj  + Cg 

D(z)  = z“"^I  + z®‘2d^_2  + . . . + zDj  + Dg 


Then  from  (7.19)  we  have 


Q(z)C(z)  + R(z)D(z)  = n(z) 
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Therefore,  if  the  system  (F,G,H)  is  N-step  reachable,  we  can  find  the 
compensator  W(.{z)  = C(z)D"^(z)  which  results  in  a closed-loop  system 
with  preassigned  dynamics.  The  realization  of  this  compensator,  which 
is  called  the  controllable  realization  ( see  FUHRMANN  [1976]  ),  is 
derived  in  the  appendix. 


CHAPTER  EIGHT 
CONCLUDING  REMARKS 


By  introducing  the  augmentation  approach  into  the  theory  of 
discrete-time  time-varying  systems  we  have  been  able  to  design  feedback 
control  systems  for  a broad  class  of  time-varying  systems  in  the  first 
part  of  this  dissertation.  Via  augmentation,  we  can  turn  an  N-step 
reachable  time-varying  system  into  a system  with  a fixed  basis,  which  is 
a very  restricted  class  of  time-varying  systems.  Using  this  fixed 
basis,  we  were  able  to  derive  a generalized  control  canonical  form  and 
the  cyclizability  result.  These  results  then  led  to  the  design  of 
feedback  control  systems  which  result  in  a closed-loop  system  matrix 
that  is  algebraically  equivalent  to  a constant  matrix  with  arbitrary 
assignable  eigenvalues.  This  assignability  result  corresponds  to  past 
work  (see  KAMEN  et  al . [1985]  and  POOLLA  [1984]),  but  our  results  yield 
controllers  having  a much  lower  order  in  general  than  the  controller 
derived  by  KAMEN  et  al . [1985]. 

In  the  second  part  of  this  dissertation  , we  developed  the 
resultant  theory  of  time-varying  systems.  From  the  discrete-time  time- 
varying  version  of  the  resultant  matrix  we  were  able  to  obtain  a test 
for  coprimeness  of  polynomials  and  polynomial  matrices  with  time-varying 
coefficients,  to  derive  a Bezout  identity,  and  to  develop  a compensator 
design  based  on  a polynomial  model.  Moreover,  we  could  relate 
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reachability  of  a system  to  the  rank  of  the  resultant,  which  resembles 
results  in  the  time-invariant  case. 

There  are  many  open  problems  which  deserve  investigation  in  the 
area  of  time-varying  systems.  We  expect  that  the  theories  developed  in 
this  dissertation  will  be  useful  in  solving  some  of  these  problems.  We 
outline  some  of  these  below  . 

3)  Continuous  time  case;  We  have  not  treated  the  continuous-time  case  in 
this  dissertation.  However,  the  augmentation  approach  may  be  applied  to 
the  continuous-time  case  easily.  For  example,  WOLOVICH  [1968]  treated 
systems  having  a fixed  lexicographic  basis,  which  belong  to  a subset  of 
the  class  of  uniformly  controllable  time-varying  systems.  But,  if  we 
use  the  augmentation  approach,  then  it  is  possible  to  design  feedback 
control  systems  for  uniformly  controllable  systems  since  the 
augmentation  results  in  a system  with  a fixed  basis.  Actually,  we  can 
design  feedback  systems  for  a broad  class  of  time-varying  systems 
including  uniformly  controllable  systems.  We  should  note  that  the 
semi  linear  transformation  which  plays  an  important  role  in  our 
development  has  a different  form  in  the  continuous-time  case.  The 
reader  may  refer  to  WOLOVICH  [1968]  and  BUCY  [1968]. 
b)  Robustness:  In  the  time-invariant  case,  even  if  the  closed-loop 
system  is  stable,  the  compensator  itself  may  not  be  stable  (see  SHAW 
[1971]).  This  may  affect  the  sensitivity  of  the  system.  A similar 
situation  may  occur  when  the  compensator  for  time-varying  systems  is 
designed  by  the  augmentation  approach.  However,  our  constructive  method 
for  the  augmentation  is  very  flexible  and  gives  infinitely  many 
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different  matrices  for  the  augmentation.  Therefore,  by  choosing  an 
appropriate  augmentation,  instability  of  the  compensator  may  be 
eliminated. 

c)  Tracking  and  disturbance  rejection:  We  have  not  treated  the  problem 
of  tracking  and  disturbance  rejection  in  this  dissertation.  However,  we 
believe  that  the  augmentation  approach  would  be  very  useful  in  studying 
this  problem.  It  is  very  interesting  that  in  the  time-invariant  case 
the  most  general  robust  controller  which  regulates  the  system  with 
disturbances  must  have  an  augmented  structure  (see  DAVISON  and 
GOLDENBERG  [1976]). 

Theory  of  the  resultant:  In  the  time-invariant  case,  the  resultant 
matrix  is  deeply  related  to  the  stability  problem  (see  JURY  [1976]). 
Therefore,  as  an  extended  project,  one  might  consider  the  relationship 
between  the  resultant  matrix  and  stability  in  the  time-varying  case.  It 
would  be  very  useful  if  we  can  extract  the  pole-zero  structure  of  time- 
varying  systems  from  the  resultant  matrix.  The  reader  may  refer  to 
KAMEN  [1986]  for  the  pole-zero  structure  of  time-varying  systems.  Other 
topics  which  one  might  consider  are  the  development  of  an  algorithm  for 
finding  common  factors  of  polynomials  and  polynomial  matrices  with  time- 
varying  matrices  and  to  extend  the  resultant  theory  of  the  discrete-time 
case  to  the  continuous-time  case. 

Systems  over  rings:  We  can  apply  the  augmentation  approach  to  the 
construction  of  dynamic  state-feedback  for  linear  systems  over 
coimnutative  rings,  if  such  commutative  rings  have  the  same  property  as 
Proposition  3.4.  For  example,  EMRE  and  KHARGONEKAR  [1982]  have  shown 
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that  coefficient-assignment  for  a finite-free  linear  system  over  a 
commutative  ring  can  be  achieved  by  dynamic  feedback  compensators.  If 
the  augmentation  approach  is  used  to  the  construction  of  compensators 
for  such  systems,  then  we  can  get  the  same  result.  Moreover,  our 
compensator  has  a much  lower  order  than  the  compensator  derived  by  EMRE 
and  KHARGONEKAR  [1982].  For  more  details,  see  KAMEN  and  HWANG  [1986], 


APPENDIX 

CONTROLLABLE  REALIZATION  OF  TIME-VARYING  SYSTEMS 

In  Chapter  Seven,  we  obtained  a compensator  that  is  described  by  a 
right  polynomial  matrix-fraction  representation.  From  this  form,  we 
will  obtain  the  controllable  realization  of  the  system  in  this 
appendix.  This  is  the  discrete-time  time-varying  version  of  FUHRMAN 
[1976]'s  standard  controllable  realization.  For  the  time-invariant 
case,  the  reader  may  refer  to  FUHRMAN  [1976]  and  BROCKETT  [1970]. 

Let  the  time-varying  transfer  function  be  given  by 

(A.l)  W(z)  = C(z)D-l(z) 

where  C(z)  e A[z]P^"^  and  D(z)  e A[z]"’X'^.  In  the  following,  we  denote 
C(z)  and  D(z)  as 

C(z)  = z'^"^C^_2  + . . . + zCi  + Cq 

(A. 2) 

D(z)  = z^I  + + . . . + zDi  + Dg 

By  a power  series  expansion  (see  Kamen  et  al . [1985]),  we  can  get 
(A. 3)  W(z)  = z"^Wj  + z“^W2  + . . . + + . . . 

where  e A. 
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(A.4)  PROPOSITION.  (KAMEN  et  al . [1985])  Let  (F,G,H)  be  a given  time- 
varying  system.  Then  the  transfer-function  matrix  W(z)  of  the  system  is 
given  by 


W(z)  = H(zl  - F)-1g  [] 

Define  the  controllable  realization  (F,G,H)  of  the  transfer 
function  W(z)  given  by  (A.l)  and  (A. 2)  as 


F 

(A.5) 


0 . • . 0 -Og 

I 

I . • . 0 

II 

0 

• • 

• 

• • 

• 

• 

I A-l. 

— 1 

o 

H = [aW^ 


Then  from  Proposition  A.4  we  must  show  that  the  realization  (F,G,H) 
given  by  (A.5)  satisfies 


W(z)  = C(z)D-^(z)  = H(zl  - F)'1g  . 


Define 


X (z)  = (zl  - F)"^G 
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and  let  X‘(z)  = [X'^^X'^  . . . X'^]  where  X^.  e A[z]^^^.  Then  we  have 


(A.6) 


zX, 


zX, 


L n J 


. 0 -Df 


-Di 


I -0. 


X > 
1— » 

1 

A 

X2 

• 

• 

• 

A 

X 

II 

J 

L n . 

L 

I 

0 


If  we 

rewrite  the  above 

equation, 

then  we 

obtai n 

A 

zXi 

A 

+ 

D.X 
0 n 

A 

= 

I 

(A. 7) 

zX^ 

• 

• 

• 

A 

+ 

DjX 

• 

• 

• 

'n-1 

• 

• 

• 

A 

zX 

n 

+ 

D ,X 
n-1  n 

X 1 

n-1 

By  combining  each  equation  in  (A. 7),  we  get 


z^X  + ,X^  + . . . + zD,X  + D.X  = I 

n n-1  n 1 n 0 n 


This  is  equivalent  to 


(A.8) 


D(z)X^  = I 


From  (A. 7)  and  (A.8),  we  have 


Ill 


X„(z)  = D-l(z) 

(A.S)  X„.i(z)  = (zl  - D„.i)D-1(z) 

• • 

• • 

• • 

X^(z)  = + . . . + Di)D"^(z) 

If  we  rewrite  (A. 9)  in  matrix  form,  then  we  obtain 


(A. 10) 


X(z)  = 


zH-li  + 2n-2g^^^  ^ 
zH-^l  + + 


. + 0]^ 
. + D2 


D-^(z) 


zl 

I 


+ D 


n-1 


We  now  consider  W(z)  = C(z)D"^(z).  By  comparing  both  sides  of 
W(z)D(z)  = C(z)  , we  have 


Cn-i  = a^Wi 

(A. 11)  Cp_2  = a'^W2  + 

• • 

• • 

• • 

Cg  = a^W|^  + (a'^  ^'■'^n-l )*^n-l  (<J^W2)D2  + (aW2)D2^ 

By  combining  (A. 10)  and  (A. 11),  we  can  get 


[aWi  a^W2  . . . a"^WjX(z)  = C(z)D-l(z) 
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From  the  definition  of  X(z)  and  (A. 5), 

H(zl  - F)‘^G  = HX(z) 

= [aWi  . . . a^WjX(z) 

Therefore,  H(zl  - F)"^G  = W(z),  and  thus  (F,G,H)  defined  by  (A. 5)  is  the 
controllable  realization  of  W(z)  = C(z)D"^(z). 

If  the  transfer  function  W(z)  is  proper,  i .e.  if  the  power  series 
expansion  of  W(z)  is  given  by 

W(z)  = Wq  + z + . . . 


then  we  have  the  controllable  realization  (F,G,H,J)  such  that  F,  G,  H 
are  defined  by  (A. 5)  and  J = Wq. 
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